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A recently  developed  filtration  model  based  on  an  effective  medium 
approximation  has  been  verified  for  Brownian  particles  in  the  presence  and  absence  of  an 
energy  barrier.  Previously,  good  agreement  between  filtration  models  and  experiments 
was  seen  only  in  the  absence  of  significant  electrostatic  repulsion,  i.e.,  when  the 
deposition  was  transport-controlled.  The  inability  to  predict  filtration  rates  when 
electrostatic  repulsion  is  significant  is  due  to  the  lack  of  an  adequate  understanding  or 
description  of  colloidal  interactions.  Therefore,  an  independent  method  of  accounting  for 
colloidal  interactions  has  been  employed  in  this  research  by  directly  measuring 
deposition  kinetics  in  a well-defined  geometry.  Computerized  video  microscopy, 
coupled  with  an  appropriate  transport  model,  was  used  to  extract  deposition  rate 
constants  between  colloidal  silica  and  glass  substrates  in  a parallel  plate  flow  cell.  In 
addition,  good  agreement  was  seen  between  the  measured  and  predicted  deposition 
behavior  in  the  parallel  plate  flow  cell  for  a variety  of  physicochemical  conditions.  So 


Xll 


the  predictions  of  the  filtration  model  based  on  an  effective  medium  approximation  for 
the  capture  of  0.5  pm-diameter  silica  spheres  in  a packed-bed  were  verified  by 
incorporating  independently  measured  deposition  rate  constants.  In  addition,  the 
filtration  of  bacteria  and  viruses  was  enhanced  via  surface  modification  of  granular  filter 
media.  However,  the  removal  rates  of  microorganisms  in  a packed-bed  could  not  be 
reliably  predicted,  thus  illustrating  the  need  for  better  understanding  of  colloidal 
interactions  in  filtration  processes. 
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CHAPTER  1 
INTRODUCTION 


1-1  General  Purpose 

The  deposition  and  dispersion  of  colloidal  particles  are  of  fundamental  interest  for 
a wide  variety  of  applications  including  drug  delivery,  printing  inks  and  municipal  water 
treatment  as  well  as  the  processing  and  manufacturing  of  microelectronic  devices, 
minerals,  cosmetics  and  biomaterials.  In  particular,  granular  and  fibrous  filtration  is  of 
practical  interest  for  removing  microorganisms  and  other  contaminants  from  municipal 
and  industrial  water  supplies  and  from  conditioned  air  as  well  as  for  improving  the 
transport  of  minerals  in  soils  for  agricultural  purposes.  Therefore,  an  adequate 
understanding  of  filtration  processes  has  significant  social,  public  health  and  economic 
implications. 

Ultimately  a methodology  for  predicting  filtration  behavior  based  on  measurable 
system  properties  is  desired.  Unfortunately,  however,  filtration  rates  for  most  natural 
systems  of  practical  interest  cannot  be  reliably  predicted  based  on  current  theories. 
Filtration  efficiencies  can  be  predicted  for  some  simple  systems  but  not  for  more 
practically  relevant  problems  such  as  the  removal  of  poliovirus  in  a packed-bed  filter. 
Therefore,  this  work  was  aimed  at  verifying  a recently  developed  filtration  model  under 
conditions,  which  have  previously  shown  poor  agreement  between  experiments  and 
models. 
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In  addition,  novel  surface  coatings  for  granular  filter  media  were  investigated  to 
enhance  the  removal  of  bacteria  and  viruses  in  a packed-bed  filter.  These  coatings  were 
often  shown  to  improve  removal  rates  of  microorganisms  by  orders  of  magnitude. 
Therefore,  advancement  in  the  technology  and  understanding  of  such  processes  can  have 
significant  societal  impacts  including  the  preservation  of  human  life. 

1-2  Accomplishments 

A filtration  based  on  an  effective  medium  approximation  was  verified  in  a 
packed-bed  filter  for  a range  of  physicochemical  conditions  using  a novel  method,  which 
eliminated  the  need  of  an  explicit  expression  for  colloidal  interactions.  By  employing  the 
interaction  force  boundary  layer  approximation,  colloidal  interaction  terms  in  the 
transport  equations  were  reduced  to  a single  parameter  in  the  boundary  condition  at  the 
surface.  This  parameter,  which  is  termed  the  deposition  rate  constant,  was  then 
independently  measured  under  well-defined  flow  conditions  in  a system  with  properties 
nearly  identical  to  those  used  in  packed-bed  filtration  experiments.  The  measured 
deposition  rate  constant  was  subsequently  used  in  conjunction  with  the  effective  medium 

approximation  to  adequately  predict  filtration  rates  for  a variety  of  physicochemical 
conditions. 

The  transport  equations,  with  the  gravity  terms  retained,  were  appropriately  scaled 
and  solved  using  a finite  difference  method  to  model  colloidal  deposition  in  a parallel 
plate  flow  cell.  The  predictions  of  this  model  were  then  compared  to  measured 
deposition  behavior  of  colloidal  silica  in  a flow  cell  observed  using  computerized  video 
microscopy.  Good  agreement  was  seen  between  predicted  and  measured  collection  rates 
of  0.5  and  1 .0  pm-diameter  silica  particles  on  glass  plates  when  deposition  was  transport- 


controlled.  For  the  cases  in  which  deposition  was  controlled  by  repulsive  electrostatic 
interactions,  collection  rates  were  much  greater  than  predicted  using  the  interaction  force 
boundary  layer  approximation  in  conjunction  with  DLVO-theory.  Deposition  rate 
constants  were  then  extracted  for  the  systems  of  interest  and  used  as  a parameter  to 
predict  filtration  rates  using  the  effective  medium  approximation.  The  system  was  also 

characterized  for  other  necessary  parameters  to  predict  filtration  rates  entirely  based  on 
existing  theories. 

In  addition,  coatings  were  applied  to  granular  filter  media  and  colloidal  silica  to 
enhance  the  inactivation  or  removal  of  microorganisms.  Metal  hydroxides  and 
quaternary  ammonium  siloxane  coatings  were  used  to  impart  a positive  charge  on  filter 
media  to  aid  in  the  capture  of  negatively  charged  bacteria  and  viruses.  Both  of  these 
coatings  were  promising  in  their  ability  to  enhance  microbe  capture,  but  the  filtration 
rates  could  not  be  reliably  predicted. 

1-3  Dissertation  Organization 

This  communication  presents  research  done  to  verify  a filtration  model  and  to 
enhance  filtration  rates  of  microorganisms.  Chapter  2 gives  a short  review  of  models  and 
experiments  found  in  the  literature  which  are  pertinent  to  this  research,  including 
colloidal  interactions  and  transport  in  a packed-bed.  Chapter  3 is  devoted  to  the 
modeling  and  experimental  efforts  for  colloidal  transport  in  a parallel  plate  flow  cell. 

Here  good  agreement  is  demonstrated  between  experimentally  measured  and 
theoretically  predicted  deposition  rates.  Chapter  4 is  dedicated  to  the  work  done  for 
verifying  the  filtration  model  based  on  the  effective  medium  approximation.  Chapter  5 
describes  all  the  critical  materials,  procedures  and  characterization  results  pertinent  to 
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deposition  experiments  in  a flow  cell  and  a packed-bed.  Finally,  Chapter  6 discusses  the 
efficacy  of  coatings  to  filter  or  inactivate  bacteria  and  viruses.  An  economic  feasibility 
study  of  aluminum  and  ferric  hydroxide  coatings  is  also  presented  here. 


CHAPTER  2 

BACKGROUND  AND  LITERATURE  REVIEW 


Predicting  colloidal  deposition  or  aggregation  requires  a description  of  the 
interaction  and  external  forces,  hydrodynamics  and  the  physical  and  chemical 
characteristics  of  the  suspension.  Much  work  has  been  done  to  understand  and 
mathematically  describe  colloidal  interactions  and  hydrodynamics,  so  a review  of 
relevant  literature  is  given  below. 

2-1  Colloidal  Interactions 

Colloidal  particles,  which  roughly  range  from  as  small  as  ten  angstrom  (10'9m)  to 
about  one  micrometer  (10‘6m)  in  diameter  (Hiemenez  and  Rajagopalan,  1997),  in 
suspensions  have  astounding  properties  due  to  the  dominance  of  interparticle  forces 
(Russel  et  al.,  1989).  Colloidal  forces  control  many  aspects  of  colloidal  suspensions 
including  their  rheological  and  optical  properties  (Russel  et  al.,  1989).  Subtle  changes  in 
physical  properties  can  lead  to  large  changes  in  colloidal  interactions  and  thus  physical 
properties.  Colloidal  interactions  have  been  divided  into  several  types  including 
electrostatic,  van  der  Waals,  hydrophobic,  hydration  and  steric  contributions. 

2-1-1  DLVO-Type  Interactions 

The  only  widely  accepted  explicit  formulation  of  these  interactions  comes  from 
the  widely  accepted  theory  of  Deijaguin-Landau-Verwey-Overbeek  (DLVO),  which 
considers  the  total  interaction  to  be  a superposition  of  electrostatic  and  van  der  Waals 
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interactions  (Hiemenez  and  Rajagopalan,  1997;  Adamson,  1990).  Independent 
developments  made  by  Derjaguin  and  Landau  (1941)  and  by  Verwey  and  Overbeek 
(1948)  form  the  basis  of  the  so-called  DLVO  theory,  which  has  since  been  modified  and 
adapted  by  many  researchers  including  Hogg  et  al.  (1966)  and  Gregory  (1981).  Several 
of  these  models  can  be  found  in  Rajagopalan  and  Kim  (1981),  Elimelech  et  al.  (1995)  and 
Li  ( 1 998).  In  addition,  a great  deal  of  experimental  work  has  been  devoted  to  measuring 
or  confirming  expressions  that  describe  electrostatic  and  van  der  Waals  interactions 
(Clapp  and  Dickinson,  2001;  Ducker  et  al.,  1994;  Meagher,  1992;  Israelachvili  and 
Adams,  1978).  For  simplicity,  only  one  model  is  used  in  this  work  to  describe  DLVO 

theory,  as  described  below.  For  constant  surface  potential,  |vy|  < 50mVand  Kap>5,  the 


electrostatic  interaction  between  a sphere  and  plate  is  given  as  (Hogg  et  al.,  1966) 


<t>e  = 2™AVoi  + X|4) 


VoiVoa 


(Voi  + ViL) 


In 


f \ + exp(-ich)^ 


1 - exp(-tch) 


+ ln(l  - exp(-2i<h)) 
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(2.1) 


Here  h is  the  minimum  separation  distance  between  the  spherical  particle  and  plate,  ap  is 
the  particle  radius,  s is  the  dielectric  constant  of  the  medium,  v|/oi  is  the  potential  of 
surface  i,  and  k is  the  reciprocal  Debye-length  defined  as 


2 e 
K = 


V n,z,2  . 


£kBT  i 


(2.2) 


The  unretarded  van  der  Waals  interaction  for  a sphere  and  a plate  is  given  by 
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where  Av  is  the  Hamaker  constant  (Hamaker,  1937).  Figure  2-1  shows  four  typical 
DLVO  interaction  profiles.  A comparison  of  models  by  Rajagopalan  and  Kim  (1981) 
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Figure  2-1.  Typical  DLVO  profiles  sphere-plate  interactions  calculated  using  a 
constant-potential  model  for  electrostatic  interactions,  which  is  represented  by  dashed 


( ) lines  and  an  unretarded  van  der  Waals  model,  which  is  represented  by 

dotted  ( ) lines.  Solid  lines  represent  the  superposition  of  electrostatic  and  van 


der  Waals  interactions.  Part  a)  represents  an  overall  attraction,  b)  shows  the  case  of  a 
large  secondary  minimum,  and  c)  represents  the  case  in  which  a barrier  to  deposition 
exists.  A is  the  Hamaker  constant,  ap  is  the  particle  radius,  c is  the  salt  concentration 
and  <))  represents  the  zeta  potential  of  the  sphere  or  plate. 
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; Figure  2-2.  The  effect  of  potassium  chloride  concentration  of  the  ratio  of 
particle  size  to  Debye  length  for  0.5  and  1 .0  pm-diameter  particles. 


suggests  that  these  expressions  are  as  likely  as  any  others  to  be  sufficient  in  describing 
deposition  phenomena.  Also,  the  requirement  that  tcap>5  holds  in  the  limit  of  large 
particles  and  small  double  layers,  and  Figure  2-2  shows  the  requirement  for  0.5  and  1.0 
pm  particles  used  in  this  experiment  is  met  at  and  above  potassium  chloride 
concentrations  of  1 0'5  M. 

2-1-2  Hydrophobic  Forces 

When  non-polar  molecules  are  introduced  into  an  aqueous  suspension,  an 
energetically  unfavorable  disruption  in  the  orientation  of  water  molecules  will  occur. 

The  system  will  tend  towards  a more  favorable  configuration,  which  can  result  in  the 
aggregation  of  the  non-polar  substrates  (Hiemenez  and  Rajagopalan,  1997;  Elimelech  et 
al.,  1995).  The  most  common  example  of  this  behavior  is  the  demixing  of  water  and  oil. 
Many  other  implications  of  hydrophobic  forces  have  been  noted  and  studied  in  recent 
years  (Israelachvili,  1992;  Farrah  et  al.,  1981;  Kokkoli,  1998).  Claesson  and  Christenson 
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(1998)  used  a surface  force  apparatus  (SFA)  to  measure  hydrophobic  forces  at  separation 
distances  up  to  80  nm  and  found  this  interaction  exceeded  van  der  Waals  interaction  by 
two  orders  of  magnitude.  Rabinovich  and  Yoon  (1994)  used  an  atomic  force  microscope 
(AFM)  to  measure  hydrophobic  forces  and  discovered  the  attractive  force  sharply 
increased  for  surfaces  with  advancing  contact  angles  greater  than  95°.  Hydrophobic 
forces  are  known  to  play  a crucial  role  in  colloidal  deposition  in  biological  systems, 
mineral  flotation,  food  processing,  etc.  and  is  one  interest  of  this  study. 

2-1-3  Hydration  Forces 

For  two  surfaces  to  make  true  contact  in  an  aqueous  solution,  all  water  adsorbed 
at  the  point  of  contacts  must  be  removed,  thereby  increasing  the  amount  of  energy 
required  to  bring  the  surfaces  into  contact.  This  additional  force  required  for  deposition 
is  known  as  hydration  forces  in  water,  and  decays  exponentially  with  a decay  length  of  a 
few  nanometers  (van  de  Ven,  1989).  While  these  forces  are  ever  present  in  colloidal 
suspension,  no  quantitative  representation  exists,  so  they  will  not  be  included  in  any 
calculations  in  this  work. 

2-1-4  Polymeric  Forces 

When  polymers  are  added  to  colloidal  suspensions,  surface  properties  and  particle 
interactions  can  be  significantly  altered  via  a variety  of  mechanisms.  Steric  stabilization 
occurs  when  adsorbed  layers  of  polymers  shield  attractive  interactions  and  also  prevent 
overlap  of  the  polymer  layers.  Bridging  can  occur  when  large  molecular  weight 
polymers  adsorb  to  more  than  one  particle  forming  an  aggregate.  Ionic  polymers  can  also 
be  used  to  induce  a certain  charge  on  a substrate  thereby  changing  colloidal  interactions. 
Polymeric  forces  are  quite  important  for  a number  of  applications  yet  little  is  know  about 
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their  origin  (Hiemenz  and  Rajagopalan.  1997;  Elimelech  et  al.,  1995),  and  they  will  not 
be  studied  in  this  work. 

2-1-5  Specific  Interactions 

Forces  such  as  electrostatic  and  van  der  Waals  are  ubiquitous  in  aqueous, 
colloidal  suspensions,  but  other  forces  such  as  receptor-ligand  interactions  are  specialized 
to  certain  molecules.  Specific  forces  are  expected  to  play  a significant  role  in  the 
adhesion  of  microbiological  particles  and  could  account  for  the  large  discrepancy 
between  theoretical  deposition  rates  and  experimental  observations  (Bell  et  al.,  1984). 
Much  effort  is  being  devoted  to  identifying,  measuring,  modeling  and  predicting  such 
specific  forces  (e.g.,  Dickinson,  1997;  Bell  et  al.,  1984),  but  they  will  not  be  explicitly 
considered  in  this  study. 

2-2  Hydrodynamics 

The  fluid  velocity  fields  required  to  solve  the  convective-diffusion  equation  (2.24) 
are  obtained  by  solving  the  equations  of  motion  for  the  appropriate  geometry.  Due  to  the 
spatial  complexity  of  a packed-bed  filter,  geometrical  approximations,  such  as  those 
discussed  in  2-3-2,  are  necessary  to  estimate  the  flow  field  around  granular  collectors. 

The  velocity  field  obtained  for  the  effective  medium  approximation  for  granular 
filtration,  as  well  as  for  flow  between  two  plates,  is  presented  in  this  section. 

2-2-1  Flow  between  Parallel  Plates 

The  two-dimensional  velocity  profile  (2.5)  for  flow  between  parallel  plates  is 
readily  obtained  by  solving  a momentum  balance  equation  with  no-slip  boundary 
conditions  (2.4a-c)  (Bird  et  al.,  1960). 

d2u  1 dP 


dy2  p dx 


(2.4a) 
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u = 0 at  y = 0 and  y = b 


(2.4b.c) 


u = 


jr  AP 
2(j.  L 


T|_  [V2 

lb J lb; 


(2.5) 


The  volumetric  flow  rate,  Q,  and  the  average  velocity,  U,  are  found  by 


b wb3AP 

Q = J wu(y)dy  = wbU  = 

o \2\iL 


(2.6) 


Also,  it  is  often  more  convenient  to  speak  in  terms  of  shear  rates  rather  than  velocity 
profiles,  so  the  shear  rate  at  the  surface  is  given  as 


(2.7) 


2-2-2  Filtration  Models 

An  adequate  model  of  the  fluid  flow  in  a packed  bed  should  describe  the 
velocity  field  around  the  collectors  and  the  pressure  drop  through  the  bed  in  a realistic 
way.  Obtaining  the  true  flow  field  in  a packed  bed  is  virtually  impossible  due  to  the 
complexity  and  irregularity  of  its  geometry,  so  geometric  models  have  been  applied  to 
approximate  the  velocity  and  pressure  fields.  The  capillary  tube  model  provides 
reasonable  prediction  for  pressure  gradients  but  not  for  the  velocity  field,  while  the 
spherical  cell  models  proposed  by  Happel  (1958)  or  Kuwabara  (1959)  provide  reasonably 
good  predictions  for  both  pressure  and  velocity  fields,  although  they  are  not  without  some 
ambiguities  (Tien,  1989;  Li,  1998).  Also,  these  cell  models  are  theoretically  limited  to  the 
case  in  which  porosity  is  above  about  0.6  (Li,  1998;  Happel  and  Brenner,  1965),  and  the 
non-uniformity  of  collector  particles  can  not  be  rigorously  taken  into  consideration  (Li 
and  Park,  2000;  Spielman,  1977).  Li  and  Park  (2000)  proposed  a filtration  model  based 
on  the  effective  medium  approximation  that  divides  the  flow  into  two  regions;  the 
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spherical  fluid  envelope  region  near  the  collector  is  described  by  creeping  flow,  and  the 
outer  or  effective  medium  region  is  described  by  Brinkman's  equation  to  account  for  the 
presence  of  other  collectors.  The  effective  medium  approximation  for  filtration  was 
found  to  be  in  better  agreement  to  experimental  permeability  measurements  and  to 
numerical  transport  calculations  than  other  models  (Li  and  Park,  2000). 

2-2-2-1  Capillary,  Hapel  and  Kuwabara  models 


a)  b) 


Figure  2-5.  Some  typical  geometrical  representations  of  granular 
packed-bed  filters:  a)  spherical  cell  model;  b)  capillary  tube  model; 
and  c)  constricted  tube  models. 


Capillary  and  cell  models,  such  as  those  developed  by  Happel  (1958)  and 
Kuwabara  (1959),  have  been  extensively  used  in  the  literature  to  describe  flow 
phenomena  in  a packed-bed  filter.  Figure  2-5  gives  a schematic  representation  of  these 
types  of  geometrical  approximations.  Capillary  tube  models  essentially  represent  a 
granular  packed-bed  as  a bundle  of  these  tubes.  Capillary  models  have  been  successfully 
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used  to  predict  pressure  drops  across  porous  media,  but  they  do  not  provide  a realistic 
representation  of  flow  fields  inside  these  porous  media  (Tien,  1989;  Li,  1998).  The  cell 
models  depicted  in  Figure  2-5,  which  are  expected  to  give  more  representative  flow  fields 
around  a granular  collector,  consist  of  a solid  collector  element  surrounded  by  a 
fluid  envelope.  Some  ambiguity  to  the  boundary  condition  exists  as  either  a stress-free  or 
vorticity-free  boundary  condition  can  be  applied  at  the  outer  edge  of  the  fluid  envelope, 
corresponding  to  the  Happel  or  Kuwabar  models,  respectively.  To  address  this 
ambiguity,  Li  and  Park  (2000)  modified  the  cell  models  for  filtration  by  employing  an 
effective  medium  approximation  as  discussed  in  2-2 -2 -2.  Only  predictions  from  the 
effective  medium  approximation  for  filtration  will  be  sought  here,  but  due  to  the 
popularity  of  these  models  they  are  worth  mentioning. 

2-2-2-2  Effective  medium  approximation 

The  effective  medium  approximation  assumes  a liquid  envelope,  whose 
dimension  is  found  from  the  porosity  of  the  bed,  and  which  immediately  surrounds  the 
spherical  collector  as  in  the  Hapel  and  Kuwabara  models.  This  envelope  is  then 
surrounded  by  a so-called  effective  medium,  which  takes  into  account  the  presence  of 
surrounding  collectors  as  shown  in  Figure  2-6.  Stokes’  equation  is  solved  inside  the  fluid 
envelope,  while  Brinkman’s  equation  is  solved  in  the  effective  medium  region,  which 
extends  to  infinity.  The  dimensionless  forms  of  these  equations,  in  which  pressure  p is 
scaled  by  pU/ac,  velocity  u is  scaled  by  U and  the  radial  coordinate  r is  scaled  by  the 
collector  radius  a c,  are  given  by  equations  (2.8)  and  (2.9): 


Vp  = V2u, 


1 <r  <R 


(2.8) 
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Figure  2-6.  Schematic  representation  of  filtration  model  based  on  the 
effective  medium  approximation. 


Vp  = V1 2u  - A.2pu , r > R (2.9) 

T T 2 2/ 

Here,  Ap  = ac  / kp , where  kp  is  the  permeability  of  the  packed-bed.  The  boundary 

conditions  imposed  are  the  no-slip  condition  at  r=l,  the  continuity  of  velocity,  stress  and 
pressure  at  r=R  and  the  velocity  far  from  the  collector  must  approach  that  of  uniform 
flow  in  the  packed-bed. 

An  analytical  solution  can  be  obtained  for  this  problem  by  using  a stream 
function,  'F,  which  is  defined  as  follows: 
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r2  sin  9 dS 


1 5T 

U Q = 

rsin&  dr 


(2.10) 


(2.11) 
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T = f(r)sin2  9. 


(2.12) 


By  applying  the  condition  that  the  velocity  must  approach  the  uniform  velocity  far  from 
the  sphere.  f(r)  in  the  outer  region  becomes 
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(2.13) 


From  the  no-slip  boundary  condition.  f(r)  is  given  in  the  inner  region  as 
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In  addition,  the  pressure  in  these  regions  can  be  expressed  as  follows: 
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(2.16) 


A,  B,  C and  D are  obtained  by  matching  the  velocity  and  stress  at  the  junction  of  these 
two  regions  (r=R)  and  are  described  by  the  following  four  equations  (2.17a-d): 
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By  matching  the  overall  porosity  to  the  size  of  the  fluid  envelope  for  a collector  size 
distribution  density  of  P(a<;),  we  determine  the  thickness  of  the  fluid  envelope  8e  by 


00 


JacP(ac)da( 


1 -8  = 
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00 


(2.18) 


f(ac  +5e)  P(ac)da, 
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Evidence  suggests  that  R should  be  taken  as  1.028(ac  + 6e)/  ac,  where  1.028  has  been 


found  empirically  (Li,  1998).  Finally,  the  permeability  can  be  written  as 


A.2  = 


a2  3(1  -e)a 
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where  D is  found  to  be 


D(R,A.p)  = --^-[6R6X3p  + 21R5A,2  - 5R4A.p  + 45R\  - 5R3A.2  + 45R3  - RA.3p  - A.2p],  (2.20) 


Q 


Q,  = 1 80R3  + 24R5A.2  - 45R4A.2  - 9A.2  + 4X\  + 1 80R4L  + 4R6X3 

(2  21' 

+ 10R3A.p  -9RA,p  -9R5A,p  +30R2A.p  -180R3A.p  ‘ 

Thus,  a complete  description  of  the  permeability  and  fluid  velocity  in  a packed-bed  of 
spherical  collectors  is  achieved  (Li  and  Park,  2000). 

2-2-3  Sedimentation 

The  effect  of  gravity  on  colloidal  transport  can  be  considered  by  first  writing  the  net 
buoyancy  force  as 

Fg  = ApVg , (2.22) 


where  Ap  represents  the  difference  in  fluid  and  particle  densities,  V the  particle  volume 
and  g is  the  acceleration  due  to  gravity.  The  settling  velocity  Us  of  the  particle  is 
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obtained  by  balancing  the  gravitational  force  with  the  drag  force  Fd  exerted  by  the  fluid 
on  a spherical  particle 


IL  = 


2g(pp-pf)ap 

9pa 


(2.23) 


Wilson  (2000)  proposed  a method  of  determining  the  bulk-settling  rate  of  particles  by 
measuring  the  absorption  at  two  different  heights  of  a colloidal  suspension  with  an  UV- 
visible  spectrophotometer  during  settling.  Basically,  the  ratio  of  the  height-interval  to  the 
time-interval  between  inflection  points  of  the  two  absorption-time  curves  gives  the 
average  settling  rate  as  depicted  in  Figure  2-4.  Wilson  (2000)  measured  the  settling  rate 
of  the  same  1 pm-diameter  silica  particles  used  in  this  work  to  be  10-4  cm/s,  while  the 
rate  estimated  from  (2.23)  was  about  5x10 3 cm/s.  The  theoretical  rates  calculated  from 

equation  (2.23)  were  used  in  obtaining  concentration  profiles  using  the  model  discussed 
in  3-2. 


Absorbance 


Figure  2-4.  Schematic  of  uv-visible  spectrophotometer  experiments  to  determine 
the  bulk  settling  rate  of  colloidal  particles  in  suspension  (Wilson  2000) 
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2-3  Colloidal  Transport 

2-3-1  General  Approaches 

Colloidal  transport  can  be  modeled  either  by  following  the  particle  trajectory 
using  a Lagrangian  viewpoint  or  by  describing  the  probability  densities  of  particles  from 
a fixed  or  Eularian  viewpoint  (Adamczyk  et  al.,  1983;  Elimelech  et  ah,  1995) . The 
Lagrangian  description  is  advantageous  in  that  the  finite-size  of  the  particles  is 
considered,  whereas  the  Eularian  description  treats  particles  as  point-masses,  which  is  the 
main  deficiency  of  continuum  theory  (Adamczyk  et  ah,  1983).  The  Eularian  perspective 
for  the  transport  of  Brownian  particles  can  offer  analytical  solutions,  which  can  be 
extended  to  heat  transfer  problems,  whereas  the  trajectory  analysis  approach  for  colloidal 
transport  requires  extensive  computation  time  due  to  the  stochastic  nature  of  Brownian 
motion.  In  this  work,  colloidal  transport  will  only  be  modeled  using  the  convective- 
diffusion  equation,  which  evolves  from  a Eularian  description,  yet  the  results  of  models 
formulated  from  the  Lagrangian  viewpoint  will  be  discussed  in  2-4. 

2-3-2  Convective-Diffusion 

A colloidal  particle  in  solution  is  subject  to  many  forces  including  Brownian 
motion,  convection,  gravity  and  colloidal  interactions,  all  of  which  are  described  by  the 
convective-diffusion  equation.  By  assuming  the  particles  will  always  tend  toward  a 
Maxwellian  distribution  of  velocities  on  an  appropriate  time-scale,  the  Fokker-Plank 
equation  was  formulated  which  describes  the  probability  distribution  of  particles 
(Chandrasekhar,  1943).  Equation  2.4  is  the  well  known  convective-diffusion  equation, 
including  the  contribution  of  external  and  surface  forces  K,  and  can  be  deduced  from  the 


Fokker-Plank  equation  for  colloidal  dispersions  in  which  particle  inertia  is  negligible 
(Adamczyk  et  al.,  1983;  Chandrasekhar,  1943). 
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dc 

h u • Vc  = V • 

dt 


DVc-D 


K 

c 

kBT 


(2.24) 


Here  c is  particle  concentration,  which  is  analogous  to  a probability  distribution,  u is 
taken  as  the  fluid  velocity,  D is  the  diffusion  coefficient,  kg  is  the  Boltzmann  constant,  T 
is  temperature,  and  t is  time.  The  Stokes  number  Nst  can  provide  information  about  the 
importance  of  inertial  impact,  even  though  no  exact  criteria  exist.  The  Stokes  number  for 
packed  beds  is  given  by 


(2.25) 


where  pp  is  the  particle  density.  Deposition  by  inertial  impaction  drastically  decreases  for 
Stokes  numbers  below  about  0.1  (Tien,  1989).  A typical  Stokes  number  for  systems  in 
this  study  is  about  0.04,  thus  justifying  the  omission  of  inertia  and  the  use  of  the 
convective-diffusion  equation  (2.24).  Difficulty  in  solving  the  convective-diffusion 
equation  arises  due  to  the  lack  of  understanding  of  colloidal  interactions  and  because  of 
the  drastically  different  length-scales  associated  with  colloidal  transport.  A method  of 
accounting  for  the  varying  length-scales,  i.e.,  the  interaction  force  boundary  layer 
approximation,  is  discussed  in  2-3-2-2,  3-1  and  3-2-1. 

2-3-2-1  Brownian  diffusion 

Brownian  motion  is  described  by  the  diffusion  coefficient,  which  is  found  by 
balancing  thermal  fluctuations  with  the  drag  exerted  on  the  particle  by  the  fluid.  Einstein 
(1920)  first  formulated  the  diffusion  coefficient  based  on  Stokes’  law  for  the  drag  on  a 
sphere  of  radius  ap: 
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Fs,  = 67t(j.apU , 


(2.26) 


where  p and  U are  the  fluid  viscosity  and  average  velocity,  respectively.  The  resulting 
equation  for  the  diffusion  coefficient  is  known  as  the  Stokes-Einstein  equation: 


kT 

D = . (2.27) 

67ipap 

However,  Stokes"  law  was  developed  for  the  drag  on  a particle  in  a fluid  media  that 
extends  to  infinity,  which  is  an  inherent  violation  in  the  deposition  process.  Brenner 
(1961)  developed  a correction  factor  to  the  Stokes-Einstein  equation  that  accounts  for  the 
effects  of  hydrodynamic  retardation  for  a sphere  moving  towards  or  away  from  a surface. 
This  correction  for  a rigid  surface  is 
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where  h is  the  separation  between  particle  and  surface  and 
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Figure  2-3  shows  the  hydrodynamic  retardation  correction  factor  as  a function  of  particle- 
wall  separation  distance  for  a 1 -micron  particle.  Employing  optical  trapping 
experiments,  Clapp  and  Dickinson  (2001)  and  Frej  and  Prieve  (1993)  have  shown  that 
Brenner’s  correction  (Eqs.  2.28  and  2.29)  for  the  position-dependent  diffusion  coefficient 
of  Brownian  particles  near  a wall  is  quite  reasonable,  especially  at  separation  distances 
large  than  the  Debye  length.  Goldman  et  al.  (1967a  and  1967b)  also  formulated  a 
correction  for  the  motion  of  a sphere  tangentially  to  the  surface,  but  this  correction  will 
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be  omitted  from  the  current  investigation  since  diffusion  tangentially  to  the  surface  has 
been  found  negligible  when  considering  the  appropriate  length-scale  in  that  direction. 


r 


Figure  2-3.  Hydrodynamic  correction  factor  for  sphere  moving  towards  a plane 
surface  (Brenner  1961). 


2-3-2- 1 Smoluchowski-Levich  solution 

The  convective-diffusion  equation  has  been  solved  for  a variety  of  geometries  and 
under  a variety  ot  simplifying  assumptions.  One  such  solution  for  purely  convective- 
diffusion  controlled  transport  of  molecular-sized  solutes  in  the  absence  of  gravity, 
hydrodynamic  retardation  and  surface  forces  is  attributed  to  Levich  (1962). 

D00V2c  = u • Vc  (2.30) 

Under  these  assumptions,  equation  2.10,  along  with  perfect  sink  and  uniform  bulk 
concentration  boundary  conditions,  Levich  ( 1 962)  found  the  surface  flux  for  flow  in  a 
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channel  to  be 


)exp[~73]d>7  bD"X 


(2.31) 


b is  the  channel  gap  and  x is  the  downstream  distance  (Levich,  1 962).  The  dimensionless 
flux  or  Sherwood  number  Sh,  defined  by  the  following  equation,  is  dependent  on 
downstream  distance  (Adamczyk  et  al.  1983;  van  de  Yen  1989): 
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(2.32) 


The  average  Sherwood  number  over  the  length  of  the  channel  L is  often  used  in  the 
literature  for  comparing  deposition  rates.  This  so-called  Smoluchowski-Levich  solution 
corresponds  to  the  diffusion-limited  case  of  particle  deposition  in  the  absence  of  gravity, 
colloidal  forces  and  hydrodynamic  interactions  and  is  equivalent  to  the  interaction  force 
boundary  layer  solution  in  the  limit  of  a large  deposition  rate  constant. 

2-3-2-2  Interaction  force  boundary  layer  approximation 

Convective-diffusion  equations  are  readily  written  to  predict  particle 


concentration  profiles  in  the  deposition  process,  but  solving  these  equations  with  gravity, 
surface  forces,  and  hydrodynamic  retardation  effects  included  becomes  quite  difficult 
even  numerically  because  of  the  drastically  different  length-scales  associated  with 
transport  and  deposition.  An  approximate  solution  can  be  obtained  by  dividing  the 
problem  into  appropriately  scaled  regions  and  matching  the  solutions  at  adjacent 
boundaries.  Using  such  an  approach,  Spielman  and  Friedlander  (1974)  introduced  the  so- 
called  interaction  force  boundary  layer  approximation  in  which  surface  force  terms 
collapse  into  an  apparent  first-order  reaction  boundary  condition  at  the  surface: 
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J = k'c , (2.33) 

ao 

k'  = D./  J{exp(<j))-l}dy  (2.34) 

0 

The  underlying  assumption  to  this  model  is  that  the  surface  force  layer  is  much  smaller 
than  the  diffusion  boundary  layer,  which  in  turn  is  much  smaller  than  the  characteristic 
length-scale  for  the  bulk  flow  domain  (e.g.  the  gap-width  between  plates  for  flow  in  a 
channel).  The  scaling  arguments  as  well  as  the  development  of  the  interaction  force 
boundary  layer  approximation  (IFBLA)  for  flow  between  parallel  plates  is  presented  in 
Chapter  3 . 

Prieve  and  Ruckenstein  (1976)  later  determined  that  the  interaction  force 
boundary  layer  approximation  originally  obtained  by  Ruckenstein  and  Prive  (1973)  and 
Spielman  and  Friedlander  (1974)  is  only  valid  in  the  presence  of  a repulsive  interaction 
between  particle  and  collector,  and  that  in  the  absence  of  an  energy  barrier,  full  transport 
equations  should  be  solved.  Using  a more  mathematically  rigorous  matched  asymptotic 
expansion  scheme,  Shapiro  et  al.  (1990)  reiterated  that  the  IFBLA  is  only  valid  in  the 
presence  of  short-ranged  repulsive  colloidal  interactions.  In  addition,  Prieve  and 
Ruckenstein  (1976)  noted  that  hydrodynamic  retardation  of  a particle  near  a wall  should 
not  be  neglected  as  it  can  reduce  the  deposition  rate  constant  by  orders  of  magnitude. 
Adamczyk  and  van  de  Ven  (1981)  also  suggest  that  the  IFBLA  is  not  appropriate  in  the 
presence  of  a significant  secondary  minimum.  Therefore,  salt  concentrations  of  colloidal 
suspension  were  kept  high  enough  to  expect  short  double  layers  but  low  enough  that  large 
secondary  minima  were  not  expected. 

The  interaction  force  boundary  layer  approximation  model  described  and  used  in 
this  work  assumes  that  a clean,  smooth,  heterogeneous  surface  is  available  for  the 
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deposition  of  non-aggregated  particles  and  that  blocking,  detachment  and  aging  effects 
are  negligible.  Arguably,  these  assumptions  hold  for  the  non-biological  experiments  in 
this  research,  in  which  particle  concentrations  in  solution  and  on  the  surfaces  are  low  and 
in  which  glass  slides  have  a root-mean-square  roughness  as  low  as  a few  nanometers 
(Truesdail,  1999).  Although  the  effects  of  blocking,  coagulation,  aging  and  detachment 
are  neglected  here,  a discussion  of  other  work  considering  these  effects  can  be  found  in 
section  2-4. 

2-3-3  Particle  Collection  by  Filtration 
2-3-3-1  Coefficient  definitions 

The  collection  of  particles  in  a packed-bed  is  usually  expressed  in  terms  of  a filter 
coefficient,  filter  filtration  efficiency,  mass  transfer  coefficient  or  Sherwood  number,  so 
definitions  of  these  terms  is  in  order.  As  usual,  the  flux  J of  particles  to  a collector  is 
given  as  the  number  of  particles  deposited  per  unit  area  per  time,  and  the  Sherwood 
number  is 


Sh  = J 


D.c0 


(2.35) 


These  terms  are  useful  in  measuring  the  local  deposition  rate,  but  the  filter  coefficient  A. 
or  efficiency  r\  better  describe  the  performance  of  the  filter.  Assuming  the  deposition 
rate  is  first  order  with  respect  to  particle  concentration,  a balance  of  particles  entering  and 
leaving  a differential  volume  of  the  column  can  be  written  as 

x t dc 

-U—  = JAsp.  (2.36) 

Here  z is  the  depth-coordinate  in  the  packed-bed,  and  the  solution  of  this  equation  yields 
the  concentration  profile  along  the  length  of  the  filter: 
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sp  m 

c(z)  = c0e  ^ Z.  (2.37) 

The  filter  coefficient,  which  is  independent  of  bed  depth,  is 

Am  km 

X = — p-— ■■ , (2.38) 

U 

Finally,  the  filtration  efficiency,  which  is  dependent  on  bed  depth  1,  is  expressed  as 

r\  = C-eff  = 1 -e~xl,  (2.39) 

cin 

where  Cjn  and  ceff  are  the  inlet  and  effluent  concentrations.  Experimentally,  the  filter 
coefficient  is  found  by  measuring  the  concentration  at  a minimum  of  two  depths  in  the 
bed,  and  equations  2.37  and  2.38  are  solved  or  fitted  to  find  the  filter  coefficient. 

2-3-3-2  Theoretical  collection  rates 

For  the  case  in  which  Brownian  diffusion,  convection  and  colloidal  interactions 
are  significant,  and  in  which  gravity,  hydrodynamic  retardation  and  interception  are 
negligible,  Spielman  and  Friedlander  (1974)  developed  a semi-analytical  solution  for  the 
deposition  rate  of  particles  onto  cylindrical  and  spherical  particles  using  the  interaction 
force  boundary  layer  approximation.  For  spherical  collectors  they  (Spielman  and 
Friedlander,  1974;  Li  and  Park,  2000)  determined  the  Sherwood  number  to  be 

Sh  = 0.635AfPel/3y,(P,),  (2.40) 

where 

Ps  =1.125As-'/3Pe~l/3Da,  (2.41) 

and  for  cylindrical  collectors,  the  Sherwood  number  was  found  to  be 

Sh  = 0.731AfPeI/3ys(ps),  (2.42) 


where 
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Ps  =1.022A;1/3Pe-,/3Da.  (2.43) 

Determining  ys  and  yc  requires  a numerical  solution  to  a Volterra  integral.  Spielman  and 
Friedlander  (1974),  therefore  tabulated  results  to  the  integral,  and  a summary  of  these 
results  is  shown  in  Table  2-1 . The  flow  parameters  As  and  Ac  depend  on  the  model  used 
to  estimate  the  flow  field  around  a collector  element,  and  they  are  summarized  for 
various  models  in  Table  2-2. 


Table  2-1.  Values  converted  from  the  Volterra  integral  solution 

by  (Spielmand  and  Friedlander,  1974) 


Ps  or  Pc 

7s 

Yc 

0 

0 

0 

0.01 

0.01404 

0.013198 

0.02 

0.02772 

0.026078 

0.05 

0.06669 

0.062857 

0.1 

0.12545 

0.118709 

0.2 

0.22410 

0.213217 

0.5 

0.42303 

0.407133 

1 

0.59840 

0.58220 

2 

0.75180 

0.73887 

5 

0.88508 

0.87783 

10 

0.93945 

0.93536 

20 

0.96895 

0.96667 

50 

0.98735 

0.98637 

100 

0.99366 

0.99307 

00 

1 

1 
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Table  2-2.  Expressions  for  the  flow  field  parameter  As  for  different  cell 
models.  Expressions  taken  from  (Li  1998;  Spielman  and  Friedlander  1974). 


Spherical  Cell 
Model 

As 

Happel 

* _ l-(l-e)5/3 

s l-(l-e)2 -X(l-s)1/3  +K(1-s)5/3 

Brinkman's 

As  = 1 + ^p 

Effective 

Medium 

Approximation 

2 |"30R(4  - R 2 ) + 1 OA.  R(-8  + 1 2R  + 5R 3 ) + 3 A2  (1  - 5R 2 1 
.A.  — v p 

5 Q,[-10R3+14R5)  + 3a3pR(1  + 5R2-10R3+4R5)  J 

where 

Q,  = 1 80R3  + 24RSX;  - 45R4A.2  - 9A.2  + 4A,3p  + 1 80R4A.p  + 4R6A,3p 
+ 1 0R3A.3p  - 9RA.p  - 9R5A.p  + 3 OR 2 A,2  - 1 80R\ 

2-4  Literature  Review  of  Deposition  Experiments 

An  extensive  amount  of  research  has  been  performed  to  understand  various 
effects  on  colloidal  deposition  in  a number  of  different  geometries  including  packed- 
beds,  rotating  discs  and  parallel  plate  flow  channels.  Studies  on  the  influence  of  physical 
(e.g.  particle  size  and  flow  rate)  and  chemical  (e.g.  salt  concentration  and  pH)  parameters 
on  particle  attachment  and  detachment  rates,  blocking  effects,  aging  phenomena,  etc.  can 
be  found  throughout  the  literature.  Unfortunately,  most  of  the  efforts  to  model  such 
effects  seem  to  be  parametric  and  non-predictive.  However,  some  of  the  results  of 
modeling  and  experimental  work  on  colloidal  deposition  is  presented  below,  and  good 
reviews  of  deposition  models  and  experiments  can  be  found  in  Elimelech  et  al.  (1995), 
Tien  (1989)  and  Adamczyk  et  al.  (1983). 
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In  general,  when  particles  and  collectors  are  oppositely  charged,  fairly  good 
agreement  has  been  repeatedly  observed  experimental  deposition  data  and  theoretical 
calculations  (Bowen  and  Epstein,  1979;  Fitzpatrick  and  Spielman,  1973;  Adamczyk  and 
van  de  Ven,  1981;  Fitzpatrick,  1972;  van  de  Ven,  1998;  Adamcz  et  al.,  1986).  This 
theoretical  and  experimental  agreement  for  the  collection  of  colloidal  particles  in  the 
presence  of  attractive  double  layers  has  been  noted  for  a number  of  geometrical  systems 
including  packed-bed  filters  (Fitzpatrick  and  Spielman,  1973;  van  de  Ven,  1998),  parallel 
plate  flow  channels  (Bowen  and  Epstein,  1979;  Adamcz  and  van  de  Ven,  1981), 
impinging  jet  stagnation-point  flows  (Adamczyk  et  al.,  1986)  and  rotating  discs  (Hull  and 
Kitchener,  1969).  However,  as  the  size  of  the  particle  being  collected  in  a packed-bed  is 
increased  beyond  approximately  a couple  of  microns  in  diameter,  the  agreement  between 
theoretical  predictions  and  experimental  data  was  seen  to  become  poor  for  even  for 
oppositely  charged  particles  and  collectors  (Elimelech,  1994;  Tien,  1989). 

Conversely,  when  collectors  and  particles  are  similarly  charged,  experimental 

collection  rates  ot  colloidal  particles  often  differ  from  theoretical  predictions  by  orders  of 

% 

magnitudes  (Marshall  and  Kitchener.  1966;  Hull  and  Kitchener,  1969;  Tobiason  and 
O’Melia,  1988;  Bowen  and  Epstein,  1979;  Sjollema  and  Busscher,  1990a).  In  such 
experiments  the  dependence  of  particle  deposition  on  pH  or  salt  concentration  could  not 
be  predicted  in  the  presence  of  double  layer  repulsion.  Typically,  the  dependence  on  pH 
was  experimentally  seen  to  be  much  more  gradual  than  predicted  by  including  DLVO- 
type  expressions  for  colloidal  interactions  in  models  for  particle  collection  in  both 
packed-bed  (Tobiason  and  O’Melia,  1988;  Tobiason,  1987;  Fitzpatrick,  1972;  Elimelech, 
1991)  and  flat  plate  systems  (Sjollema  and  Busscher,  1989  and  1990a).  Therefore,  the 
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predicted  deposition  rates  are  sometimes  orders  of  magnitudes  different  from 
experimentally  observed  deposition  rates  (Adamczyk  et  al.,  1983;  Elimelech  et  al.,  1995). 
It  should  also  be  noted  that  Chang  and  Whang  ( 1 997  and  1 998)  recently  used  the 
Kuwabara  cell  model  along  with  a trajectory  analysis  approach  to  predict  adequately 
predict  filtration  rates  of  Brownian  particles  in  the  presence  of  repulsive  double  layer 
interactions.  However,  the  results  of  these  simulations  when  deposition  was  transport- 
controlled  provided  poorer  agreement  than  previous  models. 

Nearly  every  communication  noting  poor  agreement  between  modeling  and 
experiments  seems  to  cite  heterogeneity  of  surface  charge  or  roughness  as  the  source  of 
the  discrepancy.  Some  models  have  been  developed  to  account  for  heterogeneity  of 
charge,  surface  structure  or  particle  size  distribution  (Song  et  al.,  1994;  Vaidyanathan, 
1986;  Vaidyanathan  and  Tien,  1991 ; Tien,  1989;  Czamecki  and  Warszynski,  1987),  but 
little  work  has  been  done  to  verify  heterogeneity  as  the  source  of  the  observed 
disagreement.  Typically,  these  models  show  that  heterogeneity  can  account  for 
experimental  and  discrepancies  by  assuming  surfaces  have  patch-wise  charges  or 
protrusions.  While  these  models  are  usually  parametric,  they  prove  that  the  heterogeneity 
explanation,  which  pervades  the  deposition  literature,  is  plausible.  Other  such 
explanations  include  hydrodynamic  and  electrokinetic  coupling,  surface  charge 
regulation  and  polymeric  configurations,  but  no  sufficient  model  exists  to  predict 
deposition  or  detachment  a priori  (van  de  Ven,  1998). 

Generally  when  no  barrier  to  deposition  is  present,  the  process  is  considered  to  be 
transport-controlled,  and  the  deposition  rates  can  usually  be  predicted  by  the  convective- 
diffusion  equation.  However,  when  attractive  colloidal  interactions  are  long-ranged. 
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collection  rates  in  a packed-bed  or  in  a stagnation  point-flow  system  were  experimentally 
and  theoretically  seen  to  be  significantly  higher  than  predicted  by  the  convective- 
diffusion  limit  (Chari  and  Rajagopalan.  1985;  Adamczyk  et  al.,  1989;  Elimelech,  1991). 
For  oppositely  charged  particles  and  collectors,  the  filtration  rates  were  seen  to  increase 
dramatically  as  the  salt  concentrations  were  decreased  (Elimelech,  1991;  Adamczyk  et 
al.,  1989),  which  is  opposite  to  the  trend  seen  for  similarly  charged  materials. 

Direct  observations  of  colloidal  deposition  in  well-defined  geometries  have  also 
elucidated  that  blocking  effects,  bond  aging  and  cleaning  methods  can  play  a significant 
role  in  the  deposition  rates.  Meinders  and  Busscher  (1993)  concluded  that  interaction 
forces  between  polystyrene  particles  and  glass  slides  changed  over  time  by  studying  the 
resident-time  dependence  of  desorption  in  a parallel  plate  flow  cell.  Statistical  analysis  of 
the  distribution  of  deposited  particles  has  revealed  that  a particle  on  the  collector  surface 
significantly  affects  the  deposition  of  other  particles  (Dabros  and  van  de  Ven,  1982  and 
1983;  Meinders  et  al.,  1992;  Sjollema  and  Busscher,  1990b).  The  deposited  particles  are 
reported  to  inhibit  deposition  in  regions  ranging  from  about  3 particle  diameters  away 
(Sjollema  and  Busscher,  1990b)  up  to  675  geometrical  cross  sections  (Meinders  et  al., 
1992).  Dabros  and  van  de  Ven  (1983)  also  note  that  the  type  of  acid  used  to  clean  glass 
slides  can  appreciably  affect  deposition  rates  onto  the  slides  even  though  the  measured 
zeta  potentials  were  similar  for  all  the  treatments.  Finally,  Luthi  et  al.  (1998)  concluded 
that  the  sites  available  for  adsorption  are  transient  after  observing  that  immediately  after 
desorption  of  latex  particles,  the  surface  lost  its  ability  to  adsorb  particles  for  hours. 

While  these  blocking  and  aging  effects  can  be  substantial  and  should  not  be  ignored,  the 
experiments  in  this  research  were  designed  to  minimize  such  effects.  Namely,  particle 
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concentrations  in  the  bulk  and  on  the  surface  were  maintained  low  to  minimize  blocking 
or  coagulation  effects,  and  salt  concentrations  used  were  moderate  to  minimize  the 
likelihood  of  significant  secondary  minima  (and  thus  detachment)  and  to  satisfy  the 
interaction  force  boundary  layer  requirement  of  small  double  layers. 


CHAPTER  3 

COLLOIDAL  DEPOSITION  IN  A PARALLEL  PLATE  FLOW  CELL 


3-1  Interaction  Force  Boundary  Layer  Approximation  Requirements 

Throughout  this  work,  the  interaction  force  boundary  layer  approximation 
(IFBLA),  introduced  in  2-2-3,  will  be  used  to  model  colloidal  deposition,  thus  the  validity 
of  this  approximation  is  shown  here.  The  primary  assumption  of  the  IFBLA  is  that  the 
surfaces  forces  act  over  a region  much  smaller  than  the  size  of  the  diffusion  boundary 
layer,  which  in  turn  is  much  smaller  than  some  characteristic  bulk  length-scale.  The 
appropriate  length-scale  for  the  surface  force  and  bulk  regions  are  taken  to  be  the  Debye 
length  (Eq.  2.2)  and  collector  size  for  filtration  or  gap  width  for  flow  cell  problems, 
respectively.  The  length  scale  for  the  diffusion  boundary  layer,  which  is  not  readily 
apparent,  is  found  by  matching  the  convective  and  diffusive  terms  of  the  convective- 
diffusion  equation  in  that  layer,  as  described  in  section  3-2-1  and  by  Levich  (1962). 
Experimental  conditions  must  be  chosen  to  satisfy  these  assumptions,  so  the  applicable 
range  of  operation  conditions  for  filtration  and  flow  cell  experiments  are  discussed 
below. 

First,  the  validity  of  the  scaling  assumption  for  flow  between  two  plates,  as  shown 
in  Figure  3-1  and  Equation  (3.1),  is  discussed. 

k_i  « §D  « b (3.1) 

For  systems  in  which  electrostatic  and  van  der  Waals  forces  are  predominate,  the  Debye 
length,  which  depends  mostly  on  the  type  and  concentration  of  electrolyte  present,  can  be 
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used  as  an  appropriate  length-scale  for  the  surface  force  region.  The  diffusion  boundary 
layer  length-scale,  which  is  found  by  matching  convective  and  diffusive  terms,  is 
determined  in  section  3-2-1  to  scale  as 

SD  ~ bPe~l/3,  (3.2) 


bulk 

region 


Figure  3-1 . Schematic  of  the  three  distinct  regions  in  a parallel  plate 
flow  cell. 


where  b is  the  gap-width  in  the  flow  cell,  and  the  Peclet  number  is 


Pe  = Ub/D , 


(3.3) 


where  U is  the  average 
IFBLA  requirement,  in 


fluid  velocity  and  D is  the  bulk  diffusion  coefficient.  Thus,  the 
terms  of  measurable  variables  becomes 


K 


-i 


« 


/Ubvl/3 


V J 


«1 


(3.4) 
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Also,  it  should  be  noted  that  So  is  simply  a length-scale  and  is  in  no  way  indicative  of  the 
exact  size  of  the  boundary  layer,  but  it  is  sometimes  referred  to  as  the  boundary  layer 
thickness. 

Figure  3-2  shows  that  the  ratio  of  boundary  layer  thickness  to  gap-width  is  small  over  a 
wide  range  of  flow  rates  and  particle  sizes,  as  required.  Silica  particles  of  0.5  and  1 .0  pm 
diameters  were  chosen  for  deposition  experiments  since  they  are  large  enough  to  be 
observed  microscopically  but  small  enough  to  exhibit  Brownian  motion.  Figures  3-3  and 
3-4  show  the  effects  of  potassium  chloride  concentration  and  flow  rate  on  the  ratio  of 
Debye  length  to  diffusion  boundary  layer  thickness  for  0.5  and  1 .0  pm-diameter  particles, 
respectively.  For  potassium  chloride  concentrations  above  about  10~4  M,  the  ratio 
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Figure  3-2.  Effect  of  particle  size  and  flow  rate  on  the  relative 
diffusion  boundary  layer  thickness. 
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becomes  very  small  for  both  0.5  and  1.0  pm-diameter  particles.  Salt  concentrations  of 
about  1 mM  were  used  in  most  experiments  to  satisfy  this  requirement  and  to  minimize 
changes  in  conductivity  when  adjusting  pH  with  the  addition  of  acid  or  base.  So  the 
assumptions  of  the  interaction  force  boundary  layer  assumption  for  the  deposition  of  0.5 
and  1 .0  pm  silica  particles  are  met  for  a wide  range  of  flow  rates  and  for  potassium 
chloride  concentrations  at  or  above  0. 1 mM. 

Similarly,  the  operating  conditions  in  a packed-bed  were  examined  to  meet  the 
boundary  layer  scaling  requirements.  For  flow  in  a packed-bed  with  collector  radii  ac.  the 
Peclet  number  is  given  as 

Pe  = Uac/D,  (3.5) 

and  the  diffusion  boundary  layer  length-scale  is  taken  to  be 

5d  ~ bPe~l/3  (3.6) 

Figure  3-5  shows  that  the  ratio  of  boundary  layer  thickness  to  collector  radius,  Pe'1/3,  is 
much  less  than  one  for  collector  radii  as  small  as  1 00  microns  and  for  superficial 
velocities  as  small  as  0.05  cm/s.  In  addition,  in  Figures  3-6  and  3-7,  the  Debye  length  is 
seen  to  be  much  smaller  than  the  diffusion  boundary  layer  thickness  for  potassium 
chloride  concentrations  of  1 0"4  M and  higher.  Therefore,  the  scaling  assumptions  of  the 
interaction  force  boundary  layer  approximation  are  met  for  a wide  range  of  experimental 
conditions  employed  in  this  work. 


Ratio  of  Debye  Length  to  Diffusion  Boundary;  Ratio  of  Debye  Length  to  Diffusion 

Layer  Thickness  j Boundary  Layer  Thickness 
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KCI  Concentration  (M) 

Figure  3-3.  Ratio  of  Debye  length  to  diffusion  boundary  layer  thickness 
for  0.5  |um  particles  in  a parallel  plate  flow  cell. 


KCI  Concentration  (M) 


Figure  3-4.  Ratio  of  Debye  length  to  diffusion  boundary  layer  thickness 
for  1.0  pm  particles  in  a parallel  plate  flow  cell. 
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! Figure  3-5.  Effect  of  collector  size  and  velocity  on  the  relative  diffusion 

| 

boundary  layer  thickness  for  a 0.5  pm  particle  in  a packed-bed 
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Figure  3-6.  Ratio  of  Debye  length  to  diffusion  boundary  layer  thickness 
for  a collector  radius  of  0.1  mm  in  a packed-bed 
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Figure  3-7.  Ratio  of  Debye  length  to  diffusion  boundary  layer 
thickness  for  a collector  radius  of  0.4  mm  in  a packed-bed 


3-2  Model  for  Deposition  in  a Parallel  Plate  Flow  Cell 

Since  the  interaction  force  boundary  layer  approximation  (IFBLA)  appears  valid 
for  nearly  all  cases  of  interest  in  this  work,  a solution  to  this  approximation  is  sought.  As 
shown  in  section  2-3-3-2,  a semi-analytical  solution  to  the  IFBLA  can  be  obtained  for  a 
variety  of  geometries  when  gravity  and  hydrodynamic  retardation  are  neglected. 

However,  for  large  or  dense  particles  gravity  should  not  be  disregarded,  and  the  inclusion 
of  hydrodynamic  retardation  in  the  interaction  force  boundary  layer  approximation  may 
reduce  the  calculated  deposition  rate  constant  by  orders  of  magnitude  (Prieve  and 
Ruckenstein  1976).  Therefore,  gravity  and  position-dependent  diffusion  coefficient  will 
be  considered  in  the  scaling  and  solution  to  the  deposition  problem  between  flat  plates. 
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3-2-1  Model  Development 

As  mentioned,  the  distance  over  which  surface  forces  act  is  typically  much 
smaller  than  the  bulk  and  diffusion  length-scales,  which  creates  difficulty  in  adequately 
solving  lor  particle  flux  and  concentration  profiles.  Therefore,  dividing  the  problem  into 
multiple  length  scales,  and  matching  concentrations  or  fluxes  at  adjacent  regions  is 
appropriate.  So  the  method  of  matched  asymptotic  expansion  is  employed  here  to  solve 
the  steady-state  convective-diffusion  equation  (2.24)  in  a parallel  plate  flow  cell. 

Variables  of  the  convective-diffusion  equation  are  made  non-dimensional  by 
defining  scales  appropriate  for  each  corresponding  region,  and  a summary  of  the 
dimensionless  variables  is  listed  in  Table  3-1.  In  all  regions,  the  particle  concentration  is 
scaled  by  the  bulk  concentration  c0,  and  interaction  potential  (J)  is  scaled  by  kBT,  and  the 
coordinate  in  the  direction  of  flow  is  scaled  by  the  gap-width  b.  The  natural  length-scale 
for  the  coordinate  normal  to  the  surface  is  taken  to  be  the  gap-width  b in  the  bulk  region 
and  the  Debye  length  k in  the  surface  force  region.  The  normal-coordinate  length-scale 
5q  in  the  diffusion  boundary  layer  is  found  by  matching  convective  and  diffusive  terms  in 
that  region.  The  fluid  velocity  in  each  region  is  scaled  by  the  average  velocity  and  by  the 
length-scale  in  the  particular  layer  relative  to  the  gap- width  as  shown  in  Table  3-1.  Then 
by  inspecting  the  dimensionless  convective-diffusion  equations  for  each  region,  leading- 
order  terms  are  retained,  and  terms  with  smaller  coefficients  are  eliminated  as 
demonstrated  in  the  following  sections. 
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Table  3-1 : List  of  dimensionless  variables  scaled  for  each  of  the  three  regions. 


Region  I:  Bulk 

Region  11:  Diffusion 
Boundary  Layer 

Region  III:  Surface  Force 
Interaction  Layer 
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3-2-1-1  Bulk  region 

Convection  can  be  reasoned  to  be  the  dominant  mechanism  for  particle  transport 
in  the  bulk  region  for  the  case  of  large  Peclet  numbers.  Usually,  a uniform  particle 
concentration  is  considered  at  the  inlet  of  the  flow  cell: 


c = 1 at  x = 0 


(3.7) 


The  appropriate  variable  transformation  of  the  convective-diffusion  equation  (2.24) 
yields  the  following  dimensionless  equation  for  two-dimensional,  steady-state  flow  in  the 
bulk  region: 


~5c  _ 1 
<3x  Pe 


32 


V 


cTc  cTc 
dx2  + dy: 


2~\  | f 


4- 


/ 


Pe 


27\ 


V 


dc  d</>  ~d*<p 
8y  dy  dy- 


+ 


J 


\J±dc 

U dy 


(3.8) 


Because  Pe»l,  (3.8)  can  be  reduced  to 
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Also,  since  the  velocity  u is  not  uniquely  zero,  (3.9)  implies  that  the  concentration  is 
uniform  in  the  bulk: 


c a 1 


(3.10) 


3-2-1-2  Diffusion  boundary  layer 

Similarly,  properly  scaling  (2.24)  for  the  diffusion  boundary  layer  yields 


dc 


u 


d2c  * dc  * d(f> 


* _._d2(j)*  „ dc 


dx*  dy  *2  dy  * dy  * 


+ c * — '-r-  + U 
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(3.11) 


where  U 

O 


— Pe2/3 . 
U 


Also,  since  the  interaction  force  laver  is  assumed  to  be  small 


compared  to  the  diffusion  boundary  layer,  the  interaction  potential  is  neglected  in  the 
diffusion  boundary  layer.  Thus,  (3.1 1)  is  reduced  to 


. d c * d2 c * TT  dc  * 

* = 1-  U 

dx  * dy  *2  8 dy  * 


(3.12) 


In  (3.12)  the  gravity  settling  term  Ug  may  or  may  not  be  negligible  depending  on  the 
relative  magnitue  of  the  settling  velocity  Us,  which  can  be  readily  calculated  from  (2.23), 

to  the  average  fluid  velocity  U in  the  bulk.  In  the  present  analyis,  -jj-  is  determined  to 

be  approximately  0(Pe'“  ) for  1 pm-diameter  silica  particles,  so  the  gravity  term  is 
retained. 

3-2-1-3  Surface  force  region 

Finally,  the  convective-diffusion  equation,  along  with  the  boundary  conditions  at 
the  inlet  and  at  the  surface  in  the  surface  force  region  can  be  recast  in  dimensionless  form 


as 
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At  y = 0,  c = 0 (3.13b) 

Atx  = 0,  c = 1 (3.13c) 


Since  the  interaction  force  boundary  layer  approximation  requires  that  Pe 
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and  since  — j-~0(Pe  *'  ),  (3.13a)  simplifies  to 
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J 


Equation  (3.14),  along  with  boundary  condition  (3.13b),  is  integrated  to  yield 


c = -exp(-<j))j(xf  Jexp(<j>)dy 


which  can  also  be  written  as 


c = - exp(-<t>) J (x)J  J(exp(<j>)  - l)dy  + y 
For  the  outer  region  of  the  surface  force  region,  i.e.  as  y ->  oo , 


c(y  -»oo)  = -J(x)|  |(exp((())-l)dy  +y 


(3.15a) 


(3.15b) 


(3.16) 


Now,  the  outer  region  of  the  surface  force  region  is  recognized  as  being  matched  with  the 
inner  region  of  the  diffusion  boundary  layer  such  that  (3.16)  can  be  written  as 


c * (y*  = 0)  = -J  * (x*) 


oo 

J(exp(<j>)-l)dy 


0 


(3.17) 
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3-2-1-4  Matched,  leading-order  problem 

The  resulting  equations  from  each  of  the  regions  are  now  matched  to  yield  one  set 
of  differential  equations  for  the  particle  concentration  in  the  diffusion  boundary  layer 
with  the  outer  boundary  condition  coming  from  the  bulk  region  and  the  surface  boundary 
condition  coming  from  the  surface  force  region.  Thus,  the  leading-order  problem  cast  in 
terms  of  the  diffusion  boundary  layer  variables,  in  which  the  velocity  profile  has  been 
linearized,  is 


5c*  5‘c* 

6y  * = , +U 

3x*  dy  * 

5c* 
8 5y* 

(3.18a) 

At  x*=0. 
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(3.1 8d) 
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OO 

k = 1/  J{exp(^)-l}dy  (3.18e) 

o 

An  analytical  solution  to  this  problem  is  not  available,  so  a numerical  solution  was  sought 
using  the  Crank-Nicolson  finite  difference  method  as  described  in  5-6  and  by  Hoffman 
(1989). 

3-2-2  Model  Predictions 

3-2 -2-1  Convective-diffusion  limit 

In  the  absence  of  repulsive  colloidal  interactions,  the  deposition  process  is 
transport-controlled,  i.e.  once  a particle  becomes  sufficiently  close  to  a surface  such  that 
attractive  surfaces  forces  are  significant,  the  particle  will  deposit.  Such  a case  is 


44 


mathematically  equivalent  to  a perfect  sink  boundary  condition  and  to  the  condition  of  a 
large  deposition  rate  constant  in  the  interaction  force  boundary  layer  approximation. 
Figure  3-8  shows  that  the  flux  predicted  by  the  1FBLA  matches  the  flux  predicted  by  the 
Smoluchowski-Levich  solution  (Eq.  2.31)  as  the  deposition  rate  constant  k approaches 
infinity  when  neglecting  gravity  and  hydrodynamic  retardation.  The  agreement  shown  in 
Figure  3-8  is  not  surprising,  but  it  lends  credibility  to  the  numerical  scheme  used  to  solve 
equations  (3.18a-e). 
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Figure  3-8.  Particle  surface  flux  in  the  absence  of  gravity  is  shown  as  a function  of  I 
downstream  distance  in  a parallel  plate  flow  cell  for  several  values  of  k.  For  large  ! 
k,  numerical  results  agree  well  with  the  analytical  solution  of  the  Smoluchowski-  | 
Levich  solution  (2.3 1 ).  i 


Since  the  density  of  the  silica  particles  used  in  this  work  is  about  twice  that  of  the 
suspending  electrolyte  solution,  the  gravity  term,  which  scales  to  the  square  of  particle 
radius  and  to  the  negative  one-third  power  of  fluid  velocity,  should  be  considered.  The 
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particle  flux  will  decrease  as  the  fluid  velocity  increases  when  gravity  is  significant,  but 
for  the  case  in  which  gravity  is  insignificant,  the  flow  rate  will  not  significantly  affect  the 
deposition  rate  to  a flat  plate,  assuming  the  hydrodynamic  force  does  not  overwhelm  the 
colloidal  interactions.  Figure  3-9  shows  that  the  effect  of  fluid  velocity  is  much  more 
significant  for  1 pm-diameter  silica  particles  than  for  0.5  pm  particles.  Flow  an  average 
fluid  velocity  above  about  1 cm/s,  the  deposition  rate  of  0.5  pm  silica  particles  is  not 
affected  by  the  flow  rate.  Such  effects  are  likely  to  be  observed  in  a packed-bed  filter,  so 
0.5  pm  silica  particles  will  be  used  for  filtration  experiments  to  make  use  of  the  semi- 

analytical  solution  of  Spielman  and  Friedlander  (1976)  in  section  2-3-3.2,  which  neglects 
gravity. 
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Figure  3-9.  The  influence  of  fluid  velocity  on  dimensionless  surface  Flux  in 
a parallel  plate  flow  cell  for  0.5  and  1 .0  micron  diameter  silica  particles  in  the 
transport-limited  regime 
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The  size  of  a colloidal  particle  drastically  affects  the  deposition  behavior  in  both 
the  transport-  and  interaction-controlled  regimes.  Figure  3-10  shows  the  effect  of  particle 
radius  on  the  flux  of  spherical  particles  with  a specific  gravity  of  2.0  in  a parallel  plate 
flow  cell  when  the  deposition  is  controlled  by  transport,  i.e.  when  no  barrier  to  interaction 
is  present  and  when  the  surface  is  considered  a perfect  sink.  For  small  particles  (e.g.,  a 
radius  of  50  nm),  diffusion  is  the  dominating  transport  mechanism,  whereas  gravity  is  the 
controlling  transport  mechanism  for  larger  particles  (e.g.,  a radius  of  0.5  pm).  Figure  3- 
10  shows  a minimum  in  the  flux  for  particles  with  a radius  of  about  0.2  pm.  No 


7.E-05 


i 0.E  + 00  T 1 ! 1 , 

! 0 0.1  0.2  0.3  0.4  0.5  0.6 


Particle  Radius  (microns) 

Figure  3-10.  Effect  of  particle  size  on  the  flux  of  spherical  particles  with  a 
specific  gravity  of  two  in  a parallel  plate  flow  cell. 


minimum  is  observed  for  neutrally  buoyant  particles  depositing  to  a flat  plate,  as 
illustrated  in  Figure  3-11.  In  a packed-bed,  however,  such  minima  are  seen  at  a particle 


47 


diameter  of  about  1 jam  even  for  neutrally  buoyant  particles  (Elimelech  et  al.  1995).  The 
reason  for  these  minima  is  that  as  the  particle  size  increases  above  about  1 pm,  the  inertia 
of  the  particle  becomes  significant,  thus  giving  an  increase  in  filtration  rates.  Therefore, 
error  bars  for  several  positions  in  each  of  these  figures  indicate  standard  deviations  of  the 
average  fluxes.  Since  the  particle  flux  is  dependent  on  position,  the  average  values 
shown  in  Figures  3-10  through  3-12  all  have  some  error  associated  with  them. 
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Figure  3-11.  Effect  of  particle  size  on  the  flux  of  spherical  particles  in  the 
absence  of  gravity. 


Also,  the  divergence  between  the  dimensional  and  dimensionless  flux  values  is 
worth  mentioning  at  this  point.  Because  the  diffusion  boundary  layer  length  scale 
contains  the  particle  size  and  fluid  velocity,  the  behavior  depicted  by  the  dimensionless 
flux  can  be  misleading.  For  instance,  in  Figure  3-11,  the  dimensionless  flux  is  plotted 
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against  the  particle  radius  for  spherical  particles  with  a specific  gravity  of  2.0,  similar  to 
Figure  3-10.  However,  no  minimum  is  observed  in  the  plot  of  dimensionless  flux 
because  the  scaling  factor  contains  the  particle  size.  Therefore,  one  should  be  cautious  in 
interpreting  the  results  expressed  in  dimensionless  quantities  as  phenomena  observable  in 
dimensional  values. 
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Figure  3-12.  The  dimensionless  flux  dependence  on  particle  radius  in  the 
presence  of  gravity. 
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Particle  concentration  profiles  in  a parallel  plate  flow  cell  can  be  valuable  in 
interpreting  the  influence  of  system  parameters,  such  as  gravity  or  the  deposition  rate 
constant,  on  deposition  behavior.  As  such.  Figures  3-13  and  3-16  compare  the  calculated 
concentration  profiles  obtained  for  0.5  pm-diameter  particles  in  the  presence  and  absence 
of  the  gravitational  settling  term  for  a large  and  small  deposition  rate  constant, 
respectively.  Figure  3-13  reveals  that  there  is  a depletion  of  particles  near  the  surface,  as 
compared  to  the  bulk,  when  the  deposition  is  transport-controlled,  whether  or  not  gravity 
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is  included.  For  particles  that  settle  due  to  gravity,  the  depletion  is  less  extreme;  in  other 
words,  the  gravitational  force  enhances  the  transport  of  particles  to  the  surface. 
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Figure  3-13.  Gravity  effect  on  the  concentration  profile  for  0.5  pm-diameter 
particles  in  the  absence  of  an  energy  barrier  using  a position-dependent 
diffusion  coefficient. 


3-2-2-2  Surface  interaction  controlled  deposition 

The  flux  of  colloidal  particles  is  reduced  in  the  presence  of  a significant  repulsive 
barrier.  The  deposition  rate  constant  of  the  interaction  force  boundary  layer 
approximation  is  indicative  of  the  magnitude  of  the  repulsive  barrier.  A large  deposition 
rate  constant  corresponds  to  the  lack  of  an  energy  barrier  between  particle  and  collector, 
and  a small  deposition  rate  constant  corresponds  to  a large  barrier  to  deposition.  Figures 
3-14  and  3-15  show  the  effect  of  the  deposition  rate  constant  on  the  flux  of  1 pm  and  0.5 


pm-diameter  silica  particles,  respectively. 
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Figure  3-14.  The  calculated  flux  of  1 .0  pm-diameter  silica  particles  for  a wide 
range  of  deposition  rate  constants  in  the  absence  of  hydrodynamic  retardation. 
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Figure  3-15.  The  calculated  flux  of  0.5  pm-diameter  silica  particles  for  a 
wide  range  of  deposition  rate  constants  in  the  absence  of  hydrodynamic 
retardation. 
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The  calculated  concentration  profile  for  a case  in  which  the  deposition  is 
controlled  by  repulsive  colloidal  interactions  is  depicted  in  Figure  3-16.  Here,  the 
accumulation  of  particles  near  the  surface  due  to  sedimentation  can  be  seen.  When 
gravity  is  neglected,  however,  only  a very  slight  depletion  near  the  surface  as  compared 
to  the  bulk  is  observed.  Also,  if  the  gravity  term  is  equivalent  to  the  deposition  rate 
constant,  the  concentration  profile  is  flat  and  the  flux  of  particles  towards  the  surface  is 
determined  by  sedimentation. 
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Figure  3-16.  Gravity  effect  on  the  concentration  profile  for  0.5  pm-diameter 
particles  in  the  presence  of  an  energy  barrier  using  a position-dependent 
diffusion  coefficient. 


The  predicted  flux  of  0.5  pm-diameter  silica  particles  midway  downstream  in  the 
flow  channel  is  shown  in  Figure  3-17  as  a function  of  the  deposition  rate  constant  k.  For 
deposition  rate  constants  near  ten,  the  flux  is  seen  to  be  transport-controlled,  and  further 


52 


increase  in  the  rate  constant  does  not  increase  the  flux.  These  results  are  compared  to 
experimental  deposition  rates  to  extract  deposition  rate  constants,  which  are  used  to 
predict  filtration  rates. 
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Figure  3-17.  Predicted  flux  of  0.5  mm-diameter  silica  particles  in  the 
center  of  the  parallel  plate  flow  channel. 
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3-2-2-3  Effect  of  hydrodynamic  retardation 

Because  diffusion  is  a dominant  mechanism  for  the  deposition  of  colloidal 
particles,  hydrodynamic  retardation  near  the  surface  is  expected  to  affect  the  rate  at  which 
particles  are  collected.  The  solution  for  the  position-dependent  diffusion  coefficient 
developed  by  Brenner  (1961)  predicts  a diffusion  coefficient  of  zero  at  the  surface,  but  a 
particle  needs  only  approach  close  enough  to  the  surface  such  that  colloidal  forces 
dominate.  Since  the  outer  limits  of  the  surface  force  region  corresponds  to  the  wall 
surface  in  the  model,  Brenner’s  correction  should  be  adjusted  such  that  the  surface  in  the 
model  corresponds  to  a non-zero  diffusion  coefficient  or  a finite  drag  force.  Therefore, 
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the  position  used  to  calculate  the  correction  factor  to  Stokes’  law  for  drag  should  be 
offset  from  the  position  used  in  the  model.  Figure  3-18  shows  the  effect  of  the  magnitude 
of  this  offset  on  the  flux  in  a parallel  plate  flow  cell.  As  the  size  of  the  offset  decreases, 
so  does  the  flux,  which  is  expected  since  an  offset  of  zero  corresponds  to  the  case  of  no 
diffusion  at  the  surface  and  since  a large  offset  corresponds  to  the  case  in  which 
hydrodynamic  retardation  is  disregarded.  For  purposes  of  comparing  experimental  and 
theoretical  results,  the  Debye  length  is  used  as  the  offset,  as  in  section  3-3. 
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Figure  3-18.  The  transport-limited  flux  of  1 mm-diameter  silica  particles  with  and 
without  the  inclusion  of  hydrodynamic  retardation  for  various  offset  value 
parameters. 


A comparison  of  the  calculated  concentration  profiles  for  constant  and  position- 
dependent  diffusion  coefficients  is  shown  in  Figure  3-19  for  1 pm-diameter  silica 
particles  in  the  absence  of  an  energy  barrier.  For  the  case  in  which  the  diffusion 
coefficient  is  taken  to  be  constant,  the  particle  concentration  is  depleted  near  the  surface 
and  gradually  increases  to  the  bulk  concentration  some  distance  away  from  the  surface. 
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However,  when  Brenner's  correction  to  Stokes'  law  is  considered,  the  particle 
concentration  has  a maximum  a small  distance  from  the  surfaces,  as  shown  in  Figure  3- 
19.  This  maximum  appears  because  hydrodynamic  retardation  causes  the  diffusion 
coefficient  to  become  very  small  near  the  surface,  resulting  in  a bottle-neck-type 
behavior.  In  other  words,  particles  far  from  the  surface  in  the  diffusion  boundary  layer 
are  able  to  diffuse  towards  the  surface  rather  quickly  but  are  slowed  near  the  surface  due 
to  the  increased  drag. 
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Figure  3-19.  Effect  of  hydrodynamic  retardation  on  concentration 
profiles  for  1 pm-diameter  particle  at  two  positions  in  the  flow  cell 


Because  an  analytical  solution  to  the  problem  of  particle  collection  in  a packed- 
bed  is  desired,  an  equivalent,  constant  correction  to  Stokes'  law  is  sought.  The  calculated 
flux  profiles  of  0.5  pm-diameter  silica  particles  in  a parallel  plate  flow  cell  are  compared 
for  various  correction  factors  in  Figure  3-20.  While  no  constant  correction  factor  gives 
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an  exact  profile  match  to  Brenner's  correction,  choosing  an  alpha  of  1.3  yields  an  average 
flux  that  is  within  two  decimal  places  of  the  average  calculated  using  the  position- 
dependent  correction.  Therefore,  the  constant  alpha  value  of  1 .3  is  used  to  calculate  filter 
coefficients  for  0.5  pm-diameter  particles  in  a packed-bed. 
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l 

I 

Figure  3-20.  Comparison  of  calculated  flux  for  0.5  pm-diamter  silica  particles  in 
a parallel  plate  flow  cell  for  constant  and  position  dependent  corrections  to 
| Stokes'  drag  on  a shpere. 


3-3  Experimental  Results 

3-3-1  Concentration  Effect 

The  flux  of  colloidal  silica  to  a flat  plate  was  seen  to  follow  first-order  deposition 
kinetics  with  respect  to  the  bulk  particle  concentration.  The  deposition  rate  of  1 pm  silica 
particles  in  a parallel  plate  flow  channel  was  seen  to  have  a linear  dependence  on  particle 
concentration  whether  the  rate  was  limited  by  transport  or  repulsive  colloidal  interactions. 
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Figure  3-21  shows  the  concentration  dependence  for  the  flux  of  1 pm  silica  particles  for 
two  different  concentrations  of  pH.  The  case  of  pH  2.7  represents  the  transport-limited 
regime  since  electrostatic  repulsion  is  negligible  near  the  isoelectric  point,  and  the  case  of 


pH  4.3  is  an  example  in  which  electrostatic  repulsion  is  significant.  The  linear 
dependence  of  particle  flux  on  concentration  seen  in  Figure  3-21  is  an  indication  that  the 
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deposition  kinetics  are  first-order  as  described  by  the  IFBLA  and  as  assumed  in  filtration 
models. 

3-3-2  Flow  Rate  Effect 
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Figure  3-22.  Comparison  of  experimental  data  to  flux  predicted  using  the 
interaction  force  boundary  layer  approximation. 


The  flux  of  1 pm  silica  particles  in  a 1 mM  potassium  chloride  solution  at  pH  2.7, 
which  is  close  to  the  isoelectric  point  of  silica  particle  and  glass  (see  Figure  5-13),  is 
shown  in  Figure  3-22.  As  predicted,  an  increase  in  flow  rate  was  seen  to  result  in  the 
decrease  of  the  flux  of  1 pm  silica  particles.  As  previously  mentioned,  when  electrostatic 
repulsion  is  negligible,  as  should  be  the  case  near  the  isoelectric  point  of  the  collector  or 
particles,  the  deposition  rate  is  expected  to  be  transport-controlled.  So,  the  data  in  Figure 
3-22  is  compared  to  the  convective-diffusion  limit  for  1 pm  particles  with  and  without 
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the  inclusion  of  hydrodynamic  retardation.  The  data  is  also  compared  to  the  convective- 
diffusion  limit  for  particles  with  a diameter  of  0.9  pm,  which  is  one  standard  deviation 
below  the  mean  particle  size.  Clearly,  the  experimental  data  is  in  reasonable  agreement 
with  predictions  based  on  the  interaction  force  boundary  layer  approximation  for  a 
deposition  rate  constant  of  1 0. 

3-3-3  Position-Dependent  Flux 

As  shown  in  Figure  3-14,  the  calculated  flux  of  1 pm-diameter  silica  particles  is 
dependent  on  the  position  in  the  flow  cell  when  a barrier  to  deposition  is  present  but 
independent  of  position  when  the  deposition  is  transport-controlled.  In  Figure  3-23, 
experimental  measurements  at  various  values  of  pH  qualitatively  verify  these  behaviors. 
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Figure  3-23.  Comparison  of  experimentaly  and  calculated  position- 
dependent  behavior  of  1 pm-diameter  silica  particles  in  a flow  cell. 
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At  pH  2.7,  which  is  near  the  isoelectric  point  of  the  silicates,  the  flux  is  seen  to  be 
relatively  independent  of  position.  At  pH  4.5  and  6,  however,  the  particle  flux  shows  a 
strong  dependence  on  position.  Theoretically,  this  behavior  provides  a good  means  of 

extracting  deposition  rate  constants,  but  the  scatter  in  the  data  is  often  too  large  to  reliably 
do  so. 

As  shown  in  Figure  3-14,  the  flux  is  not  expected  to  vary  with  downstream 
position  in  the  flow  cell  when  deposition  is  transport  controlled.  In  general,  such 
behavior  was  observed  experimentally,  but  a fairly  large  scatter  in  the  data  was  seen.  For 
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Figure  3-24.  Position  dependence  of  deviation  from  the  average  flux  for  the 
deposition  of  1 pm-diameter  silica  particles  at  pH  2. 7-2. 8. 
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1 pm-diameter  silica  particles  at  pH  2. 7-2. 8,  the  deviation  of  individual  flux 
measurements  from  the  average  of  these  measurements  is  shown  in  Figure  3-24.  Here, 
the  deviation  or  residual  value  is  generally  seen  to  be  smaller  downstream  in  the  flow  cell 
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than  near  the  entrance.  Also,  the  residuals  are  not  seen  to  be  either  positively  or 
negatively  biased. 

3-3-4  Dependence  on  pH 

As  the  pH  is  increased  from  the  isoelectric  points  of  the  colloidal  particles  and 
glass  substrates,  repulsive  electrostatic  barrier  between  the  particle  and  substrate  becomes 
increasingly  significant.  Therefore,  the  flux  of  silica  particles  is  expected  to  decrease  as 
pH  is  increased  from  the  transport-limited  case  at  low  pH  near  the  isoelectric  point. 


x 

3 


(0 

(A 

aj 

c 

o 

'35 

c 

d) 

E 


2.5 


1.5 

1 


£ 0.5 
0 


DLVO  prediction 
■ Experimental 


6 

pH 


8 


9 

10 


Figure  3-25.  The  flux  of  1 mm-diameter  silica  particle  averaged  for  several 
positions  in  a parallel  plate  flow  cell.  The  experimental  dependence  on  pH  is 
much  more  gradual  than  that  predicted  by  DLVO  theory  based  on  a fit  of  the 


Figures  3-25  and  3-26  show  that  the  experimentally  measured  particle  flux  in  a parallel 
plate  flow  channel  gradually  decreases  as  the  pH  is  increased  for  1 .0  and  0.5  pm- 
diameter  silica  particles,  respectively.  The  gradual  decrease  is  in  contrast  to  the  drastic 
reduction  predicted  by  DLV O theory  based  on  a fit  of  the  measured  zeta  potentials,  as 
seen  in  Figure  3-25.  Figure  3-26  shows  the  effect  of  pH  on  the  flux  0.5  pm-diameter 
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Figure  3-26.  The  effect  of  pH  on  the  flux  of  0.5  mm-diameter 
silica  particles  in  a parallel  plate  flow  cell. 


silica  particles  in  the  center  of  the  flow  channel.  As  mentioned  in  2-4,  this  is  in  keeping 
with  the  findings  of  other  researchers  for  a wide  variety  of  geometries.  By  comparing  the 
results  of  Figure  3-26  to  predicted  deposition  rates  in  Figure  3-17,  deposition  rate 
constants  were  extracted,  and  then  these  extracted  rate  constants  were  used  to  predict 
filtration  rates  with  the  effective  medium  approximation  for  comparison  to  experimental 
filter  coefficients  in  section  4-3. 


CHAPTER  4 

VERIFICATION  OF  FILTRATION  MODEL 


The  filtration  of  colloidal  silica  has  been  carried  out  in  columns  packed  with 
spherical  glass  beads  under  a variety  of  physicochemical  conditions.  The  filtration 
coefficients  obtained  from  these  experiments  were  compared  to  a filtration  model  based 
on  the  effective  medium  approximation,  as  described  by  Li  and  Park  (2000)  and  in 
section  2-3-3.  Good  agreement  was  seen  between  predictions  and  experimental 
observations  over  a wide  range  of  pH  when  independently  measured  deposition  rate 
constants  were  applied  to  the  effective  medium  approximation,  as  opposed  to  using 
deposition  rate  constants  predicted  by  DLVO  theory. 

4-1  Predictions  from  Effective  Medium  Approximation 
Filtration  coefficients  for  various  flow  rates,  collector  sizes  and  deposition  rate 
constants  were  calculated  using  the  effective  medium  approximation  model,  Happel’s 
model  and  Brinkman's  model.  Figure  4-1  shows  the  effect  of  collector  size  and 
deposition  rate  constant  on  the  filter  coefficient  for  1 pm-diameter  particles  for  all  three 
models.  As  expected,  smaller  collector  sizes  correspond  higher  filtration  rates,  and  a 
large  deposition  rate  corresponds  to  a high  filtration  rate.  Figure  4-1  also  illustrates  that 
once  the  deposition  rate  constant  exceeds  a certain  value,  filtration  can  not  be  enhanced 
with  further  increase  in  the  rate  constant.  The  predicted  effect  of  superficial  velocity  on 
filtration  rates  is  shown  in  Figure  4-2  for  multiple  values  of  the  deposition  rate  constant. 
Figure  4-3  shows  the  effect  of  pH  on  filter  coefficients  for  Happel’s,  Brinkman’s  and  the 
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effective  medium  approximation  models  using  DLVO  theory  to  calculate  deposition  rate 
constants  based  on  measured  zeta  potentials,  which  are  shown  in  Figure  5-13.  The 
predicted  filtration  rates  in  Figure  4-3  are  seen  to  have  a very  sharp  dependence  on  pH  for 
all  the  filtration  models,  and  the  models  can  not  be  distinguished  above  about  pH  3.9.  A 
more  detailed  analysis  of  the  effect  of  physicochemical  parameters  on  predictions  based 
on  the  effective  medium  approximation  can  be  found  in  Li  and  Park  (2000)  and  Li 
(1998). 


Deposition  rate  constant  k 

Figure  4- 1 . Effect  of  deposition  rate  constant  on  filter  coefficient  of  0.5  pm- 
diameter  particles  collected  by  both  0.15  and  0.3  mm-radius  spheres  in  a 
packed  bed  at  a superficial  velocity  of  1 mm/s. 


{ 


64 


3.00E+00 


MB 


E 2.50E+00  i 


£ 2.00E+00 

<D 

■ MM 

o 

■ MM 

St 

<D 
O 

o 


0) 


t 1.50E+00 


u=1mm/s,  Happel 
>u=1mm/s,  EMA 
" u=1mm/s,  Brinkman 
" — u=5mm/s,  Happel 
— u=5mm/s,  EMA 
* * u=5mm/s,  Brinkman 

u=1cm/s.  Happel 

u=1cm/s,  EMA 

u=1cm/s,  Brinkman 


1 .00E+00 


if  5.00E-01 


0.00E+00 


0.000001  0.0001 


0.01 


1 


100 


10000 


Deposition  rate  constant 

| 

i 

Figure  4-2.  The  effect  of  superficial  velocity  and  deposition  rate  constant  on 
I filter  coefficient  for  0.5  mm-diameter  particles  in  a packed-bed  of  0.6  mm- 
| diameter  spherical  collectors. 
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Figure  4-3.  Effect  of  pH  on  filter  coefficient  of  0.5  pm  particles  in  a 
packed-bed  of  0.6  mm-diamter  glass  spheres  using  DLVO  calculations 
based  on  measured  zeta  potentials. 
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4-2  Experimental  Results 

The  filtration  rates  of  0.5  and  1 .0  pm-diameter  silica  particles  have  been 
measured  for  a wide  range  of  pH,  as  well  as  for  multiple  flow  rates  and  potassium 
chloride  concentrations.  To  obtain  filtration  coefficients,  particle  solutions  were  passed 
through  a granular  packed-bed  as  described  in  section  5-1,  and  the  particle  concentrations 
were  measured  at  various  column  depths.  The  filter  coefficient  was  found  by  fitting  the 
particle  concentration  over  the  length  of  the  column  to  an  exponential  decay,  as  shown  in 
Figure  4-4.  Figure  4-5  shows  the  effect  of  fluid  velocity  on  filter  coefficient  for  1 pm- 
diameter  silica  particles  in  a packed-bed  of  0.3  mm-diameter  glass  spheres.  As  expected 
from  Figure  4-2,  an  increase  in  superficial  velocity  results  in  a decrease  in  filter 
coefficient. 

The  influence  of  salt  concentration  and  pH  on  the  filtration  of  colloidal  silica  in  a 
packed-bed  is  shown  in  Figures  4-6,  4-7  and  4-8.  Figures  4-6  and  4-7  show  the  effect  of 
potassium  chloride  and  pH,  respectively,  on  the  filter  coefficient  of  1 pm-diameter  silica 
by  0.3  mm-diameter  collectors.  An  increase  in  potassium  chloride  concentration  results 
in  a compression  of  the  double  layer,  thus  decreasing  the  range  of  electrostatic  repulsion 
thereby  causing  van  der  Waals  attraction  to  become  dominant.  Therefore,  as  shown  in 
Figure  4-6,  the  filter  coefficient  increases  as  the  salt  concentration  is  increased.  An 
increase  in  pH  is  seen  (Figure  4-7)  to  cause  a decrease  in  filter  coefficient,  since 
electrostatic  repulsion  dominates  as  the  zeta  potential  of  the  particles  and  collectors 
becomes  more  negative.  Similarly,  the  effect  of  pH  on  the  filtration  of  0.5  pm-diameter 
silica  in  a packed-bed  of  0.6  mm-diameter  collectors  is  seen  in  Figure  4-8.  The  increase 
of  filtration  coefficients  with  the  increase  in  electrolyte  concentration  and  with  the 
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decrease  in  pH  is  expected,  but  the  effects  are  much  more  gradual  than  predicted  by 
DLVO  theory.  Since  DLVO  theory  does  not  yield  deposition  rate  constants  that 
adequately  predict  filtration  rates  in  the  presence  of  electrostatic  barriers,  an  independent 
method  was  devised  to  verify  the  effective  medium  approximation  in  the  presence  of 
energy  barriers  by  independently  measuring  the  deposition  rate  constants. 
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Figure  4-4.  Silica  concentrations  relative  to  inlet  concentration  after  passing 
i through  columns  packed  with  0.6  mm  glass  beads.  Filter  coefficients  are 
determined  from  an  exponential  fit  of  the  data. 


1 

j |_ 

CM  O 

c 

i .2 

8 4- 

o 

j £ 
0) 

j 

6 j 

o 
! O 

C 

j 

4 i 

o 

■ mm 

cS 

1 

j 

2 1 

I = 

oi 

0 


♦ 

♦ 


♦ 


0.1  0.2  0.3  0.4  0.5 


♦ 

— 1 - j 

0.6  0.7 


Superficial  velocity  (cm/s) 

Figure  4-5.  The  effect  of  flow  rate  on  the  filter  coefficient  for  the  collection 
of  1 .0  mm-diameter  silica  particles  in  a packed-bed  of  0.3  mm-diameter 
glass  spheres  at  pH  4 with  a potassium  chloride  concentration  of  0.1M. 
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Figure  4-6.  The  effect  of  potassium  chloride  concentration  on  the  filter 
coefficient  for  the  collection  of  1 .0  pm-diameter  silica  particles  in  a packed-bed 
of  0.3  mm-diameter  glass  spheres  at  pH  7 and  at  a superficial  velocity  of  4 mm's. 
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Figure  4-7.  The  effect  of  pH  on  the  filter  coefficient  for  the  collection  of  1.0 

pm-diameter  silica  particles  in  a packed-bed  of  0.3  mm-diameter  glass  spheres 
at  a superficial  velocity  of  0.7  mm/s. 
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Figure  4-8.  Experimental  filtration  coefficients  for  0.5  pm- diameter  silica  particles 
in  a packed-bed  of  0.6  mm-diamter  glass  spheres  at  a superficial  velocity  of  0. 1 
mm/s. 
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4-3  Verification  of  the  Effective  Medium  Approximation 

By  using  video  microscopy  to  directly  measure  the  deposition  rates  of  0.5  pm- 
diameter  silica  particles  to  a flat  glass  slide,  deposition  rate  constants  were  extracted  with 
the  aid  ot  the  model  discussed  in  section  3-1 . The  deposition  rate  constants  are  expected 
to  describe  the  interaction  between  the  particle  and  surface  and  not  depend  on  system 
parameters,  such  as  geometry  or  flow  rate.  Therefore,  the  deposition  rate  constant 
measured  between  silica  and  glass  in  a flow  cell  are  extended  to  predict  filtration  rates 
between  colloidal  silica  and  spherical  glass  collectors  in  a packed-bed. 

The  experimental  deposition  data  from  Figure  3-26,  along  with  the  theoretical 
predictions  from  Figure  3-17,  produce  deposition  rate  constants  for  0.5  pm-diameter 
silica  particles  and  glass  at  several  pH  values.  Figure  4-9  compares  experimental 
filtration  rates  to  the  predictions  from  the  effective  medium  approximation  model  using 
both  the  measured  and  DLV O-predicted  deposition  rate  constants.  Clearly,  using  this 
independent  measure  of  the  deposition  rate  constant  gives  much  better  agreement 
between  predicted  and  measured  filter  coefficients  for  silica  in  a packed-bed  of  glass 
spheres.  Also,  Figure  4-10  shows  the  ratio  of  experimental  to  theoretical  filtration 
coefficients  for  0.5  and  1.0  pm-diameter  silica  particles  for  several  flow  rates  when 
deposition  is  expected  to  be  transport-controlled.  The  effective  medium  approximation 
predicts  the  expected  rates  within  about  a factor  of  two,  which  is  about  as  good  as  should 
be  expected  (Fitzpatrick  1972).  Therefore,  the  filtration  model  based  on  the  effective 
medium  approximation  is  seen  to  be  valid  for  the  collection  of  colloidal  particles  in  a 
packed-bed  filter. 
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Figure  4-9.  Comparison  of  experimental  filtration  rates  to  those  predicted  by  DLVO- 
theory  and  by  independently  measured  deposition  rate  constants  for  0.5  pm-diameter 
silica  particles  in  a packed-bed  of  0.6  mm-diamter  glass  spheres. 
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Figure  4-10.  The  ratio  of  experimental  and  theoretical  filtration  coefficients 
based  on  the  effective  medium  approximation  for  conditions  that  are  expected  to 
be  transport-controlled. 


CHAPTER  5 

MATERIALS,  METHODS  AND  CHARACTERIZATION 


Video  microscopy  has  been  used  to  measure  particle  deposition  in  a parallel  plate 
flow  cell  to  verily  a filtration  model  and  to  investigate  colloidal  interactions.  Colloidal 
particles,  granular  media  and  glass  slides  have  been  treated,  characterized  and  employed 
in  deposition  experiments  in  this  research.  Thus  this  chapter  describes  experimental  and 
numerical  procedures  used  and  communicates  results  from  characterization  experiments. 
Specifically,  coating  techniques,  counting  procedures,  electrokinetic  measurements, 
filtration  and  flow  cell  experiments  will  be  discussed. 

5-1  Filtration  Experiments 

Packed-bed  filtration  experiments  were  conducted  using  granular  media 
consisting  of  spherical  glass  beads,  Ottawa  sand  or  activated  carbon  to  remove  a variety 
of  colloidal  particles  including  silica,  bacteria  and  viruses.  In  all  cases,  filters  were 
charged  from  the  top  with  aqueous  colloidal  dispersionss  under  laminar  flow  conditions 
using  peristaltic  pumps.  Particulate  feed  solutions  were  diluted  to  concentrations  ranging 
from  parts  per  million  to  parts  per  thousand  to  minimize  any  particle-particle  interactions 
and  to  avoid  filter  saturation,  i.e.  to  operate  in  initial  clean-bed  stages. 

5-1-1  Filter  Configuration 

The  filter  flow  system  was  configured  for  the  convenience  of  filling,  rinsing,  and 
sampling  that  are  necessary  for  measuring  filter  coefficients  for  various  physical  and 
chemical  conditions.  Figure  5-1  shows  a schematic  of  a single-column  filtration  system. 
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which  allows  the  fluid  to  be  pumped  from  the  bottom  or  from  the  top  using  a single 
peristaltic  pump.  Filtration  columns  consisted  of  polymethylmethacrylate  (PMMA) 
tubing  (0.75  in  ID/1 .0  in  OD)  and  3 in  cylinder  end  pieces  as  shown  in  Figure  5-2.  End 
pieces  consisted  of  a disc  and  a cylinder  bolted  together  and  secured  to  the  column  by 
sandwiching  an  o-ring  between  them,  as  shown  in  Figure  5-3.  Fine  steel  screens  (about 
1 00  mesh)  are  secured  at  the  end  of  the  tube  to  prevent  granular  media  from  escaping 
from  the  column.  To  provide  a means  for  attaching  tubing,  a thread  (1/8  in.  NPT)  was 
tapped  into  the  end  piece.  All  solutions  for  filtration  experiments  were  pumped  through 
Masterflex  Tygon®  food  grade  tubing  of  1/8  in.  diameter  with  Masterflex  peristaltic 
pumps,  both  of  which  were  obtained  from  Cole-Parmer. 

5-1-2  Filtration  Procedure 

Columns  were  gradually  filled  with  granular  media  using  a tap-and-fill 
methodology.  Prior  to  adding  granular  media,  the  filtration  system  was  completely 
assembled  with  the  exception  of  the  top  end  piece,  and  a few  inches  of  water  was  charged 
into  the  bottom  of  the  column.  Then,  approximately  1 0-20  ml  of  granular  media  was 
added  to  the  top  of  the  column,  and  the  column  was  tapped  until  the  media  no  longer 
settled  or  shifted;  this  procedure  was  repeated  until  the  column  was  filled.  The  top  end 
piece  was  then  secured  to  the  column,  and  water  was  pumped  upward  until  all  the  air  was 
removed  from  the  column  and  tubing.  Often,  the  flow  direction  had  to  be  reversed  to 
facilitate  the  removal  of  trapped  air  in  the  column  and  tubing.  Once  the  column  was 
filled  and  air  removed  from  the  lines,  flow  was  directed  downward  through  the  column, 
and  the  flow  rate  was  measured  and  adjusted  as  necessary.  The  source  was  then  changed 
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Figure  5-1 . Schematic  of  packed-bed  filtration  system.  Dotted  lines  represent  path  of 

fluid  in  the  upward  flow  system. 


Figure  5-2.  Column  used  for  packed-bed  filtration  experiments. 


Figure  5-3.  Unassembled  column  used  for  packed-bed  filtration  experiments 
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from  purified  water  to  an  aqueous  hydrochloric  acid  solution  of  pH  2. 6-3.0  to  minimize 
pH  fluctuations  as  recommended  by  Elimelech  et  al.  (1995).  The  packed-bed  was  acid 
washed  in  this  manner  for  about  1 5 minutes,  after  which  the  source  was  switched  to  a 
particle-free  solution  adjusted  to  the  desired  pH  and  salt  concentration.  To  promote 
equilibrium  in  the  column,  this  particle-free  water  was  passed  through  the  column  for 
about  15  minutes  prior  to  charging  the  solution  with  particles. 

To  obtain  the  desired  particle  concentration,  typically  a few  milliliters  of  a 
concentrated  solution  would  be  added  to  about  2 liters  of  the  electrolyte  adjusted  water. 
After  particles  were  introduced  into  the  system,  some  time  was  required  for  them  to  pass 
through  the  column  and  tubing.  So,  initially  the  effluent  particle  concentration  was  zero 
and  increased  to  an  equilibrium  value  until  saturation  of  the  column  occurred.  A plot  of 
effluent  concentration  versus  time,  as  seen  in  Figure  5-4,  revealed  that  breakthrough  can 
be  confidently  expected  after  two  residence  times.  After  this  time,  the  samples  were 
taken  from  the  entrance  and  exit  of  the  column,  and  the  concentration  was  measured 
using  a UV-vis  spectrophotometer  as  described  in  section  5-3-2 -2.  Then  the  filter 
effluent  was  collected  in  a 2L  Erlenmeyer  flask  for  use  as  the  influent  in  a second  pass. 

A typical  experimental  run  consisted  of  multiple  concentration  measurements  from  the 
entrance  and  exit  of  the  column  for  two  passes,  each  of  which  was  operated  for  about  45 
minutes  to  an  hour. 

5-1-3  Granular  Media 

A variety  of  granular  media  was  employed  in  packed-bed  filtration  experiments  to 
remove  colloidal  particles  from  aqueous  streams.  Soda-lime  glass  spheres  of  various 
sizes  were  obtained  from  Mo-Sci  Corp,  Ottawa  sand  was  obtained  from  Fisher  Scientific, 
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Figure  5-4.  Typical  breakthrough  curve  for  filtration  experiments. 


and  granular  activated  carbon  was  obtained  from  U.S.  Filter  brand  cartridge  filters. 
Before  filtration,  granular  media  was  sieved  to  the  desired  size  range,  which  was  25x30 
mesh  (0.6-0.7mm  diameter)  for  most  experiments.  Once  sieved,  granules  were  soaked 
overnight  in  Chromerge  acid  solution  and  rinsed  with  deionized  water  in  an  aspirator 

until  the  pH  of  the  effluent  reached  that  of  the  water.  The  granules  were  then  dried  in  an 
oven. 

5-2  Parallel  Plate  Flow  Cell  Experiments 

Computerized  video  microscopy  was  used  to  measure  colloidal  deposition 
kinetics  in  a parallel  plate  flow  cell.  The  deposition  of  0.5  and  1.0  pm  diameter  silica 
spheres  in  a parallel  plate  flow  cell  was  directly  measured  using  video  microscopy.  The 
video  microscopy  technique,  which  was  originally  used  by  Dabros  and  van  de  Ven 
(1998),  provided  a means  of  directly  monitoring  particle  deposition  under  well-defined 
laminar  flow  conditions. 
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Figure  5-5.  Parallel  plate  flow  cell  apparatus. 


5-2-1  Equipment 

The  parallel  plate  flow  cell  obtained  from  CytoDyne  consists  of  a Plexiglas  base 
channeled  for  fluid  flow,  a glass  slide  onto  which  particles  were  deposited  and  a silicone 
gasket  to  separate  the  two  as  shown  in  Figure  5-5.  The  silicone  gasket  was  placed  over  a 
groove  through  which  a vacuum  is  pulled  to  hold  the  base  and  slide  together.  The 
dimensions  of  the  flow  domain  in  the  cell  are  8 mm  (w),  37  mm  (1)  and  254  pm  (h).  The 
glass  side  of  the  flow  cell  was  placed  downward  on  the  motorized  stage  of  an  inverted 
microscope,  and  images  of  deposited  particles  were  taken  at  various  positions  and  times 
and  stored  on  a computer.  A CCD  camera  was  used  to  view  images  from  the  inverted 
Nikon  microscope  with  a motorized  stage  that  could  be  controlled  with  a joystick  or  with 
software.  Optimas  software  package  was  used  for  controlling  the  stage,  capturing  images 
and  counting  particles.  Ruta  (1998)  gives  a detailed  explanation  of  the  automation  and 
image  analysis  program  Optimas  and  includes  the  codes  written  for  automation  and 
analysis. 
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5-2-2  Experimental  Procedure 

Immediately  preceding  deposition  experiments,  the  flow  cell,  spacer  gasket  and 
glass  slide  were  all  rinsed  with  ethanol  to  minimize  contamination.  Two  syringes  were 
used  as  teeds  to  the  flow  cell.  One  syringe  contained  the  desired  particle  solution,  and 
the  other  contained  a particle-free  solution  of  the  same  pH  and  salt  concentration.  After 
the  flow  cell  was  assembled,  a vacuum  pump  was  used  to  secure  the  slide  to  the  base,  and 
the  particle-free  solution  was  pumped  through  the  cell  at  full  speed  until  all  air  bubbles 
were  evacuated.  This  solution  was  then  pumped  through  the  flow  cell  at  the  desired  flow 
rate  for  about  1 5 minutes  to  ensure  equilibrium  and  to  inspect  the  system  through  the 
microscope.  The  strength  of  the  vacuum  was  adjusted  to  minimize  water  leakage  into  the 
pump  trap  while  keeping  an  adequate  seal.  Finally,  the  source  was  switched  to  the 
desired  colloidal  suspension  for  deposition  experiments. 

Microscope  objective  lenses  with  magnifications  of  40x  and  lOOx  were  used  to 
capture  images  for  1.0  and  0.5  pm  diameter  silica  particles,  respectively.  Since  small 
fields  result  in  higher  statistical  uncertainty,  objective  lenses  of  the  smallest 
magnification  possible  were  used.  Sample  images  taken  at  subsequent  times  of  1 pm 
silica  spheres  using  a 40x  objective  lens  are  shown  in  Figure  5-6.  Initially,  the  lower 
right-hand  comer  of  the  cell  was  marked  as  the  origin  of  the  coordinate  system.  Several 
fields  were  then  chosen  and  focused  for  image  analysis.  The  software  package  was  used 
to  recursively  scan  the  flow  cell  to  capture  images  at  several  locations  over  time.  Particle 
counts  were  seen  to  increase  linearly  in  time  as  seen  in  Figure  5-7,  and  saturation  was  not 
reached  during  the  course  of  the  experiments.  The  deposition  rate  was  determined 
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from  a least  squares  fit  of  the  slope,  and  the  particle  fluxes  were  obtained  by  dividing  the 
slope  by  the  field  area. 


Figure  5-6.  Examples  of  images  taken  for  one-micron  silica  particles  using  a 40x 
objective  lens. 
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Figure  5-7.  Typical  particle  counts  taken  at  several  locations  over  time 


79 

5-3  Preparation  and  Counting  of  Particle  Suspensions 

Dilute  colloidal  dispersions  were  prepared  for  use  in  filtration  and  deposition 
experiments  to  verify  a filtration  model.  Accurate  particle  counts  of  these  dispersions 
were  necessary  to  quantify  experimental  results  and  compare  them  to  models.  The 
procedures  used  for  preparing  and  counting  inorganic  dispersions  are  described  below, 
whereas  the  preparation  and  counting  of  biological  suspensions  are  located  in  6-1 . 

5-3-1  Dispersion  Preparation 

Colloidal  dispersions  were  prepared  at  low  concentrations  to  minimize  any 
particle-particle  interactions  or  agglomerate  formations.  Suspensions  of  0.5  and  1.0  pm- 
diameter  silica  particles  were  prepared  at  concentrations  ranging  from  0.0001-0.01 
percent  by  weight,  depending  on  the  physicochemical  conditions  of  the  dispersion  and 
the  nature  of  the  experiments.  Stock  solutions  of  about  1 % by  weight  were  made  to 
facilitate  the  preparation  of  desired  suspension  concentrations.  The  stock  solutions  were 
prepared  by  adding  a known  amount  of  dry  particles  to  membrane-filtered  water. 
Nanopure  and  Milli-Q  membrane  filtration  systems  yielded  highly  pure  water  with 
resistances  of  1 8.2-1 8.3  MQ.  The  water  was  passed  through  a 0.2  pm  filter  to  further 
remove  biological  contamination.  Prior  to  preparing  desired  suspensions,  the  stock 
solution  was  sonicated  for  at  least  four  hours  to  disperse  the  particles  and  to  break 
agglomerates.  About  one  to  three  liters  of  membrane  filtered  water  was  adjusted  to  the 
desired  salt  concentrations  and  pH.  Potassium  chloride  was  used  to  adjust  salt 
concentrations  to  the  desired  level,  which  was  usually  1 mM.  The  solution  pH  was 
adjusted  using  hydrochloric  acid  or  sodium  hydroxide  as  necessary.  Most  of  the 
electrolyte-adjusted  water  was  used  to  prepare  dilute  suspensions  from  concentrated  stock 
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solutions,  and  the  remainder  was  used  for  rinsing  and  bringing  systems  to  equilibrium. 

Suspensions  were  sonicated  for  about  one  hour,  just  prior  to  experiments  to  ensure 
dispersion. 

5-3-2  Counting  Techniques 

Reliable  particle  counting  techniques  were  needed  to  measure  deposition  rates. 

An  indirect  method  was  used  for  filtration  experiments,  whereas  a direct  counting  method 
was  used  for  parallel  plate  flow  cell  experiments.  Direct  counting  techniques  are 
advantageous  in  that  images  of  the  surface  can  be  stored  for  detailed  analysis  of 
coagulation,  blocking,  and  transport  behavior.  Indirect  counting  techniques  are  faster  and 
can  be  implemented  for  on-line  analysis. 

5-3-2-1  Direct  particle  counts 

For  parallel-plate  flow  cell  measurements,  absolute  particle  concentrations  were 
needed  since  fluxes  are  dependent  on  the  number  of  particles  in  the  bulk.  Particle  counts 
were  thus  made  directly  under  a microscope  in  a Petroff-Hauser  counting  chamber.  The 
counting  chamber  consisted  of  a grooved  glass  plate  topped  with  a microscope  coverslip. 
Since  the  spacing  in  the  groove  was  precisely  known,  two-dimensional  particle  counts 
were  converted  to  suspension  concentrations. 

5-3-2-2  Indirect  particle  counts 

In  a packed-bed  filter,  particle  concentrations  were  measured  using  an  UV-Visible 
spectrophotometer  since  number  counts  were  not  required,  and  direct  counting  techniques 
were  too  tedious  for  the  number  of  samples  to  be  counted.  A calibration  curve,  as  in 
Figure  5-8,  was  constructed  to  determine  the  weight  concentration  of  particles  from 
absorbance  measurements.  Figure  5-8  shows  a linear  relationship  between  particle 
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concentration  and  light  absorbance  at  low  particle  concentrations  (c<0.05),  as  predicted 
by  the  Beer-Lambert  law  (Alberty  and  Silbey,  1995).  When  operating  in  this  linear 
regime,  relative  particle  concentrations  were  directly  obtained  from  absorbance  ratios. 
All  experiments  were  conducted  below  about  0.03  weight  percent  to  avoid  calculating 
absolute  concentrations,  which  are  not  needed  for  filtration  experiments. 
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Figure  5-8.  Absorbance  of  350nm  wavelength  light  for  the  desired  concentration  range: 
a)  one  micron  diameter  silica  solutions  show  linear  dependence  on  absorption  for 
concentrations  up  to  about  0.04%;  b)absorbance  for  half  micron  diameter  silica  solutions, 
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5-4  Surface  Modification 

V arious  surface  properties  of  particles  and  collectors  can  be  tailored  by  physical 
or  chemical  modification  of  the  surfaces,  thereby  altering  colloidal  interactions  and 
deposition  behavior.  In  this  research,  granular  media  have  been  coated  by  physical 
adsorption  of  metal  hydroxides  and  chemical  binding  of  quaternary  ammonium  siloxanes 
to  alter  surface  properties  of  the  media.  Both  types  of  coating  were  used  to  impart  a 
positive  charge  on  negatively  charged  filter  media  to  enhance  the  capture  of  bacteria, 
viruses  and  silica,  which  are  all  negatively  charged  at  natural  pH. 

5-4-1  Metal  Hydroxide  Coatings 

Metal  hydroxides  precipitated  onto  granular  media  have  been  successfully  used  to 
enhance  the  removal  of  a wide  range  of  microorganisms  and  metals  (Lukasik  1998;  Chen, 
1998;  Scott  2000).  Lukasik  (1998)  investigated  the  efficacy  of  numerous  metal  mixtures 
and  precipitation  methods  to  enhance  the  capture  of  microorganisms.  A combination  of 
aluminum  and  ferric  hydroxide  proved  to  be  the  most  effective  and  the  procedure  for 
precipitating  this  mixture  was  used  in  this  work  (Lukasik  1 998).  An  aqueous  mixture  of 
0.25  M ferric  chloride  and  0.5  M aluminum  chloride  were  added  to  the  granular  media 
and  heated  to  60°C  until  dry.  The  media  was  mixed  twice  with  3N  ammonium  hydroxide 
and  dried  at  60°C  to  ensure  a sufficient  coating. 

5-4-2  Quaternary  Ammonium  Siloxane  Coatings 

Chemical  modification  of  silicates  (e.g.  glass  and  sand)  via  siloxane  bridging  can 
be  advantageous  to  precipitation  because  the  siloxanes  can  be  readily  functionalized  and, 
in  principle,  are  more  stable  than  physically  adsorbed  materials.  In  this  research,  two 
types  of  quaternary  ammonium  siloxanes  have  been  used  to  impart  hydrophobicity  and  a 
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positive  charge  on  the  surface  of  glass,  sand  or  silica.  Both  coatings  were  taken  from 
commercial  applications  and  are  thus  expected  to  be  safe  for  use  as  a coating  for  filtration 
media.  The  first  coating  procedure  is  referred  to  as  the  Suhara  coating  and  has  been 
employed  for  cosmetics,  and  the  second  procedure  yields  a coating  referred  to  as  AEM 
5700  and  has  been  used  commercially  as  an  antimicrobial  agent. 

5-4-2-1  Suhara  coating 

The  Suhara  coating  procedure  is  completed  in  a sequence  of  three  steps  as 
described  by  Suhara  (1995a),  and  the  reagent  in  the  final  reaction  determines  the  length 
of  the  hydrocarbon  chain  on  the  quaternary  ammonium  group.  The  flexibility  in 
hydrocarbon  chain-length  is  desired  as  the  chain-length  determines,  in  part,  the  degree  of 
the  coating  hydrophobicity.  The  chemical  structures  of  the  reagents  and  products  for 
each  reaction  is  shown  in  Figure  5-9,  where  the  subscript  m represents  the  number  of 
carbon  atoms.  The  procedure  involves  three  distinct  reactions:  (1)  polymerization  of 
2,4,6,8-tetramethylcyclotetrasiloxane  (TMCTS)  to  form  a coating  of  polymethylsiloxane 
(PMS),  (2)  hydrosilylation  with  chloromethylstyrene  (CMS),  and  (3)  introduction  of  a 
quaternary  ammonium  group  by  subsequent  reaction  with  N,N-dimethyl-n-alklyamine. 
The  details  of  each  reaction  are  described  below. 

The  first  reaction  (polymerization)  was  carried  out  either  in  a dry  or  an  aqueous 
environment.  For  the  dry  reaction  method,  about  1%  wt.  of  TMCTS  was  placed  in  a 
vessel  with  the  glass  beads  and/or  slides  and  cooked  in  an  oven  for  about  20  hours  with 
the  container  tightly  sealed  and  about  1 0 hours  with  the  container  opened.  The  aqueous 
method  was  carried  out  in  a 95%  solution  of  ethanol  adjusted  to  a pH  of  4.5-5.5  with 
acetic  acid.  Then  TMCTS  was  added  to  a concentration  of  about  2%.  The  solution  was 
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Figure  5-9.  Three  step  Suhara  coating  procedure;  a)  reaction  1 -Polymerization  of 
TMCTS;  b)  reaction  2-hydrosylation  of  PMS-Si02;  c)  reaction  3-addition  of  quaternary 
ammonium  group  to  CMS-Si02 
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stirred  for  about  5 minutes  to  allow  hydrolysis  to  take  place,  and  then  the  beads  or  slides 
were  added.  The  reaction  was  carried  out  for  about  2 or  3 minutes  and  subsequently 
rinsed  with  ethanol  and  left  to  dry  at  room  temperature  overnight. 

For  the  second  reaction  (hydrosilylation),  toluene  was  used  as  the  solvent  and 
tributylyamine  and  hydrogen  hexachloroplatinate  (IV)  hydrate  (platinic  acid)  as  catalysts. 
The  solvent,  catalysts,  chloromethylstyrene,  and  PMS  coated  glass  beads  and/or  slides 
were  refluxed  in  a sealed  container  at  1 1 0°C  for  about  4-6  hours  and  subsequently 
washed  with  toluene.  The  PMS-CMS  coated  beads  were  then  left  to  dry  overnight. 

For  the  final  reaction,  glass  beads  coated  with  PMS  and  chloromethylstyrene  were 
refluxed  with  N,N-dimethyl-«-alklyamine  and  acetonitrile  at  65  °C  for  about  6 hours. 

The  coated  beads  were  then  rinsed  with  acetonitrile  and  dried  in  vacuum  at  1 1 0°C  for 
about  1 2 hours.  N,N-dimethylbutylamine,  N,N-dimethyloctylamine  and  N,N- 
dimethyldodecylamine  were  used  as  reagents  to  yield  hydrocarbon  chain-lengths  of  4,  8 
and  1 2,  respectively. 

The  Suhara  coating  procedure  described  is  lengthy  and  tedious,  and  spectroscopic 
analysis  is  required  to  ensure  each  reaction  occurs  as  expected.  Therefore,  a simpler 
procedure  is  desired  to  increase  the  likelihood  of  chemically  binding  a quaternary 
ammonium  alkylamine  to  the  surface. 

S-4-2-2  AEM  5700  coating 

Granular  and  colloidal  particles,  as  well  as  glass  slides,  were  coated  with 
octadecyldimethyl(3-trimethoxysilyl-propyl)ammonium  chloride,  which  is  available 
commercially  from  Dow  Chemical  in  an  alcohol  solution  known  as  AEM  5700.  The 
quaternary  ammonium  siloxane  coated  onto  Ottawa  sand,  glass  and  silica  is  depicted  in 
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Figure  5-10  and  is  referred  to  as  AEM  5700  coating  in  this  communication.  A lwt% 
aqueous  solution  of  octadecyldimethyl(3-trimethoxysilyl-propyl)ammonium  chloride  was 
stirred  at  room  temperature  for  about  ten  minutes  to  allow  for  hydrolysis  and 
condensation  prior  to  adding  media.  About  200  ml  of  granular  media  (or  about  1 gram  of 
colloidal  particles)  were  then  added  and  stirred  in  the  coating  solution  for  ten  minutes. 

The  coated  media  was  removed,  rinsed  in  a vacuum  filtration  system  to  remove  any 
residual  siloxane  and  dried  overnight.  Scanning  electron  micrographs  reveal  a non- 
uniform  siloxane  coating  as  shown  in  Figure  5-11. 
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Figure  5-10.  Structure  of  3-(trimethoxysilyl)propyldimethyloctadecyl  ammonium 

chloride  (AEM  5700). 


5-5  Characterization  Techniques 

Systems  employed  in  filtration  and  deposition  experiments  were  characterized  to 
obtain  necessary  parameters  for  theoretical  predictions.  Zeta  potential  was  measured  for 
colloidal  and  granular  particles,  hydrophobicity  was  quantified  for  granular  and  flat 
substrates,  the  particle  and  collector  size  was  measured,  and  the  solution  pH  and 
conductivity  were  measured. 
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Figure  5-11.  SEM  micrograph  of  1 micron-diameter  bare  silica  (left)  and  siloxane  coated 
silica  (right). 


5-5-1  Electrokinetic  Measurements 

Electrokinetic  measurements  performed  to  determine  the  zeta  potential  of 
granular  and  colloidal  particles.  Electrophoresis  was  used  to  measure  the  zeta  potential 
with  the  commercially  available  Brookhaven  Zeta  Plus,  and  a streaming  potential  device 
was  constructed  to  measure  the  zeta  potential  of  granular  media.  Measured  zeta 
potentials  were  used  to  evaluate  coating  procedures  and  to  calculate  colloidal  interactions 
based  on  DLVO  theory. 

A streaming  potential  device  was  constructed  using  filtration  columns  adapted  to 
include  a silver  mesh  electrode  as  shown  in  Figure  5-12.  The  electrodes  were  made  of  40 
mesh  silver  gauze  welded  to  1 8 gauge  silver  wire.  The  electrodes  were  then  electroplated 
in  a 0.1  M hydrochloric  acid  solution  by  applying  a 10  mA  current  for  about  30-60 
minutes  using  a copper  cathode.  The  mesh  electrodes  also  served  to  retain  the  granular 


media  inside  the  column. 
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The  zeta  potential  was  measured  by  passing  an  electrolyte  solution  through  the 
column  at  various  pressure  drops  and  measuring  the  response  in  voltage  drop.  Potassium 
chloride  solutions  were  pumped  through  columns  of  about  20-40  cm  in  length  using  a 
peristaltic  pump,  and  the  pressure  drop  was  measured  using  either  a mercury  manometer 
or  an  inverted  water  manometer.  The  voltage  drops  were  measured  using  a handheld 
multimeter  attached  to  electrodes  at  both  ends  of  the  column.  The  voltage  drop  was 
plotted  against  pressure  drop,  and  the  slope  was  used  to  calculate  the  zeta  potential  as 

C = (AV  n K)/(AP  80  8w),  (5- 1 ) 

where  rj,  K and  sw  are  the  viscosity,  conductivity  and  dielectric  constant  of  the  electrolyte 

solution  and  s0  is  the  permitivity  of  vacuum  (Hunter,  1981;  Van  Wagenen  and  Andrade, 
1989;  Truesdail  et  al.,  1998b). 


Figure  5-12.  In-house  built  streaming  potential  device:  a)  column;  b)electrode. 
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5-5-1-1  Unmodified  media 

The  zeta  potential  of  Ottawa  sand  in  a 1 mM  potassium  chloride  solution  was 
measured  over  a pH  range  of  2.6-10.6.  Figure  5-13  illustrates  that  the  isoelectric  point  of 
the  sand  occurs  at  about  pH  2.6,  and  the  zeta  potential  of  the  sand  becomes  more  negative 
as  the  pH  is  increased.  Figure  5- 14a  shows  the  effect  of  pH  on  the  measured  zeta 
potential  of  glass  spheres  in  a 1 mM  KC1  solution.  The  isoelectric  point  of  glass  beads 
was  seen  to  be  about  2.7.  A quadratic  equation  was  fitted  to  the  data  for  interpolation  of 
zeta  potential  values  over  a wide  range  of  pH.  Zeta  potentials  of  1 pm  silica  particles 
were  measured  by  electrophoresis  and  are  shown  in  Figure  5- 14b.  Zeta  potential  for 
silica  particles  was  seen  to  be  similar  to  those  obtained  using  streaming  potential 
measurements  of  glass. 


Figure  5-13.  Zeta  potential  of  Ottawa  sand  that  has  been  soaked  overnight 
in  hydrochloric  acid. 
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Figure  5-14.  Electrokinetic  measurement  of  uncoated  silicates:  a)streaming  potential  of 
soda-lime  glass  beads;  b)electrophoresis  of  colloidal  Stober  silica. 
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5-5-1-2  Modified  media 

Streaming  potential  experiments  were  conducted  to  determine  the  zeta  potential  of 
Ottawa  sand  coated  with  metal  hydroxides  and  AEM  5700.  Figure  5-15  shows  the 
measured  zeta  potential  of  aluminum  and  iron  hydroxide  coated  sand  between  pH  values 
of  4 and  10.  The  isoelectric  point  was  determined  to  be  about  pH  8.  The  zeta  potential  of 
AEM  5700  coated  sand  was  seen  to  be  about  60  mV  immediately  after  coating  but 
quickly  dropped  and  became  negative  after  about  50  hours  of  washing  as  depicted  in 
Figure  5-16.  The  reduction  in  zeta  potential  of  quaternary  ammonium  siloxane  coated 
media  during  washing  has  also  been  observed  by  Truesdail  (1999). 
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| Figure  5-15.  Zeta  potential  of  metal  hydroxide  coated  sand. 
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Figure  5-16.  Zeta  potential  of  AEM  5700  coated  Ottawa  sand  as 
it  is  washed  with  deionized  water.  The  coated  media  loses  its 
positive  charge  quite  rapidly. 


5-5-2  Hydrophobicity 

The  hydrophobicity  of  glass  slides  and  granular  media  was  quantified  by  contact 
angle  and  capillary  uptake  measurements,  respectively.  Contact  angles  of  water  droplets 
on  glass  slides  were  determined  using  a contact  angle  goniometer.  Large  contact  angles 
corresponded  to  a high  degree  of  hydrophobicity.  To  quantify  the  hydrophobicity  of 
granular  media,  capillary  uptake  experiments  were  performed  in  which  granular  media 
was  packed  in  a glass  tube  and  placed  in  a water  bath,  as  depicted  in  Figure  5-17.  The 
level  of  water  inside  the  tube,  relative  to  the  water  bath,  was  an  indication  of  the 
hydrophobicity.  Uptake  heights  were  large  for  hydrophilic  media  and  small  for 
hydrophobic  media.  These  two  techniques  appeared  to  give  a good  indication  of  the 
degree  of  hydrophobicity  of  glass  beads  and  glass  slides. 
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Figure  5-17.  Capillary  uptake  experiments  used  to  quantify  hydrophobicity. 


Glass  beads  and  slides  were  coated  with  quaternary  ammonium  siloxanes  of 
varying  hydrocarbon  chain  lengths  using  the  Suhara  procedure  described  in  5-4-2- 1. 
Figure  5-18  shows  contact  angle  and  capillary  uptake  measurements  for  coatings  of 
varying  hydrocarbon  chain  length.  Both  types  of  measurements  indicate  that  longer 
chain  lengths  are  more  hydrophobic  than  shorter  chain  lengths. 

5-5-3  Coating  Stability 

The  longevity  or  stability  of  surface  coatings  is  of  practical  and  fundamental 
interests  and  directly  affects  the  efficacy  of  the  coatings  to  be  employed  in  filtration 
processes.  The  zeta  potential  was  measured  over  time  as  the  media  was  washed  with 
water.  Hydrophobicity  of  the  media  was  also  measured  for  the  hydrophobic  siloxane 
coatings.  Experiments  have  indicated  that  metal  hydroxide  coating  lifetimes  are  on  the 
order  of  months  (Chen  et  al.  1998)  and  that  siloxane  coatings  quickly  lose  their  positive 
charge  but  not  hydrophobicity. 
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Figure  5-18.  The  relative  hydrophobicity  of  glass  substrates  treated  with  quaternary 
ammonium  siloxanes  increases  with  hydrocarbon  chain  length,  a)  contact  angle 
measurements  of  water  on  glass  slides,  and  b)  water  uptake  heights  for  glass  beads  in  a 
1 cm  glass  tube. 
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Figure  5-19.  Stability  of  AEM  5700  coated  Ottawa  sand  that  is  washed  with  water  for 

one  month  in  a column  or  soaked  in  water  for  two  months:  a)  capillary  uptake 

experiments  show  that  the  coated  sand  retains  its  hydrophobic  behavior,  while  b) 

streaming  potential  measurements  indicate  the  positive  charge  of  the  coated  media  is  not 
retained. 


The  stability  of  metal  hydroxide  coatings  was  investigated  by  Chen  (1998)  and 
Chen  et  al.  ( 1 998)  and  found  to  be  about  4 months.  During  this  time  the  zeta  potential 
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decreased  from  positive  to  negative,  which  correlated  to  deterioration  in  performance  of 

the  filter  media.  Electrokinetic  measurements  indicate  that  the  zeta  potential  of  AEM 

5700  coated  sand  reverses  much  quicker  than  for  metal  hydroxide  coated  sand  as  shown 

in  Figure  5-16.  However,  the  coated  media  remained  significantly  more  hydrophobic 

than  uncoated  media  for  the  duration  of  the  experiments,  as  seen  in  Figure  5-19,  which 

compares  the  capillary  uptake  to  zeta  potential  for  AEM  5700  coated  sand  that  has  been 
washed. 

5-6  Numerical  Methods 

Since  an  analytical  solution  to  equations  (3-18a-e)  is  not  readily  available,  a 
numerical  solution  was  sought  to  reduce  the  differential  equations  to  algebraic  ones.  An 
implicit  finite  difference  scheme  known  as  the  Crank-Nicolson  method  was  used  to 
simplify  the  calculation  procedure  without  introducing  large  error  terms.  The  calculation 
domain  typically  contained  1 0,000  by  1 0,000  calculation  nodes.  By  using  the  Crank- 
Nicolson  scheme,  calculations  were  made  for  all  positions  of  y at  a single  x position, 
based  on  the  results  from  the  previous  x position,  i.e.  calculations  were  swept  across  the 

x-direction.  The  development  of  this  convenient  calculation  scheme  is  discussed  below, 
and  more  details  can  be  found  in  Hoffman  (1989). 

The  Crank-Nicolson  calculation  scheme  employs  central  differencing 
approximations  in  both  independent  variables.  The  x-direction  was  divided  into  intervals 
of  (Ax/2),  and  the  y-direction  was  divided  into  intervals  of  (Ay).  Nodes  in  the  y-direction 
are  denoted  by  i,  which  ranges  from  1 to  I,  and  nodes  in  the  x-direction  are  denoted  by  n. 
which  ranges  from  1 to  N.  Central  differencing  equations  were  written  for  both  variables 
at  nodal  positions  of  i and  n + !4  resulting  in  error  terms  on  the  order  of  (Ax)2  and  (Ay)2. 
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The  difference  equations  for  derivatives  with  respect  to  y were  approximated  by 
averaging  about  n and  n+1  to  eliminate  the  n + */2  node  from  the  equations.  The  resulting 
system  of  linear  algebraic  equations  readily  conformed  into  systems  of  tridiagonal 

equations,  which  were  decomposed  into  upper  and  lower  diagonal  matrices  to  further 
simplify  calculations. 


CHAPTER  6 

FILTRATION  AND  INACTIVATION  OF  MICROORGANISMS 

Removing  microbial  and  chemical  contaminants  from  municipal  water  supplies  is 
a major  public  health  concern  all  over  the  world.  A 1998  USA  Today  special  report  on 
the  dangers  of  drinking  water  cites  a Centers  for  Disease  Control  and  Prevention  (CDC) 
record  stating  that  408,000  people  became  sick  between  1993  and  1996  due  to  52 
confirmed  waterborne  illness  outbreaks  (Eisler,  1998).  One  of  the  outbreaks  in  1993  was 
linked  to  1 1 1 deaths  and  403,000  illnesses  in  the  Milwaukee  area  (Eisler,  1998).  The 
article  also  cites  a joint  report  from  the  Medical  College  of  Wisconsin  and  the  U.S. 
Environmental  Protection  Agency  (EP A)  that  suggests  7. 1 million  people  experience 
nausea  or  diarrhea  and  that  as  many  as  1,200  people  die  each  year  in  the  U.S.  from 
contaminated  water  (Eisler,  1 998).  While  accurately  confirming  the  statistics  of  illnesses 
and  fatalities  due  to  inadequately  treated  municipal  water  supplies  is  virtually  impossible, 
the  need  for  improved  processes  to  provide  clean,  safe  drinking  water  is  clear. 

Water  treatment  plants  use  a variety  of  methods,  including  chemical  treatment  and 
deep-bed  filtration,  to  treat  water  (J.  M.  Montgomery,  Inc.,  1985).  A number  of 
investigations  have  been  performed  to  improve  the  ability  of  packed-bed  filters  to  remove 
microorganisms  (Brown  et  al.,  1974;  Lukasik  et  al.,  1996  and  1999;  Chen,  1998; 
Truesdail,  1998a;  Scott,  2000;  Shaw  et  ah,  2000).  Granular  filter  media  has  been  coated 
with  peroxides  (Scott,  2000),  metal  oxides  (Lukasik  et  ah,  1996  and  1999;  Chen,  1998; 
Truesdail,  1998a;  Scott,  2000;  Shaw  et  ah,  2000)  and  organosilicates  (Fuchs,  1987)  to 
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successfully  enhance  the  adsorption  or  inactivation  of  a variety  of  microorganisms 
including  protozoa,  bacteria  and  viruses.  For  instance,  Lukasik  et  al.  (1996  and  1999) 
showed  that  virus  removal  rates  were  greater  than  99%  in  columns  packed  with  ferric  and 
aluminum  hydroxide  coated  sand,  as  compared  to  less  than  30%  removals  in  columns 
packed  with  uncoated  sand.  Here,  the  ability  to  enhance  the  removal  of  bacteria  and 
viruses  in  column  filtration  and  batch  inactivation  experiments  was  investigated  for 

media  coated  with  ferric  and  aluminum  hydroxides  and  with  quaternary  ammonium 
siloxanes. 

6-1  Materials  Used  for  Microbe  Adsorption  Studies 

Ferric  and  aluminum  hydroxides  were  precipitated  onto  granular  Ottawa  sand  and 
grade  1 diatomaceous  earth  (Sigma  Chemicals  Co.,  St.  Louis,  MO),  and 
octadecyldimethyl(3-trimethoxysilyl-propyl)ammonium  chloride  (AEM  5700)  was 
coated  onto  granular  Ottawa  sand  and  colloidal  silica  to  investigate  the  ability  of  different 
coatings  and  media  to  adsorb  or  inactivate  E.  coli,  MS-2  and  PRD-1.  Before  coating, 
Ottawa  sand  was  sieved  to  about  0.6-0. 7 mm,  which  is  a typical  size  used  in  large-scale 
municipal  and  industrial  water  filters  (J.  M.  Montgomery,  Inc.,  1985).  Granular  activated 
carbon  was  taken  directly  from  a CGAC  Series  U.S.  Filter  brand  cartridge,  sieved  to 
remove  fines  and  packed  in  columns  for  direct  comparison  to  coated  and  uncoated  sand- 
packed  columns.  Activated  carbon  was  chosen  as  a control  filter  media  because  it  is  a 
standard  component  of  point-of-use  filtration  systems  (e.g.,  household  faucet  filters  and 
water  pitchers),  which  are  becoming  increasingly  popular  due  to  concern  over 
microbiological  qualities  of  water  delivered  by  municipal  utilities.  Phages  were  assayed 
as  plaque  forming  units  (PFU)  using  agar  overlay  procedure:  MS-2  (ATCC  15597-B1) 
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was  assayed  with  Escherichia  coli  C-300  (ATCC  15597)  as  a host,  and  PRD-1  was 
assayed  using  Salmonella  typhimurium  (ATC  19585)  (Scott,  2000).  Escherichia  coli 

(ATC  13706)  was  assayed  as  colony  forming  units  (CFU)  using  a spread  plate  technique 
on  MacConkey  Agar  (Scott,  2000). 

6-2  Batch  Inactivation  of  Bacteria  and  Viruses 

Batch  studies  were  carried  out  in  which  solutions  of  coated  and  uncoated  colloidal 
and  granular  particles  were  charged  with  microorganisms,  agitated  and  assayed  for 
microbial  activity.  Ottawa  sand  that  was  unmodified,  precipitated  with  ferric  and 
aluminum  hydroxides  (as  described  in  5-4-1)  or  coated  with  AEM  5700  (5-4-2.2). 
spherical  glass  beads  that  were  unmodified,  coated  with  AEM  5700  or  coated  with  PMS- 
CMS-nn-dimethyoctylamine  (5-4-2. 1),  diatomaceous  earth  coated  with  ferric  and 
aluminum  hydroxides,  and  1 .5  pm-diameter  silica  that  was  unmodified  or  coated  with 
AEM  5700  were  each  individually  studied  for  their  ability  adsorb  and/or  inactivate  MS-2 
or  E.  coli.  The  particles  were  added  to  10  ml  of  deionized  water  in  a 50  ml  centrifuge 
tube,  and,  for  the  case  of  colloidal  silica,  were  gently  sonicated  in  a water  bath  for  30 
minutes  to  facilitate  dispersion.  The  10  ml  of  buffer,  consisting  of  0.02  M imidizole  and 
0.02  M glycine,  were  seeded  with  10  pi  of  concentrated  MS-2  or  E.  coli,  and  the  solution 
was  assayed  to  obtain  original  microbe  counts.  Then  the  microbe  solution  was  added  to 
the  centrifuge  tubes  containing  the  granular  or  colloidal  particles.  The  centrifuge  tubes 
were  mounted  on  a rocking  shaker  and  agitated  for  45-60  minutes.  After  agitation,  silica 
solutions  were  settled  by  low-speed  (1000  rpm)  centrifugation  for  5 minutes,  and 
granular  media  was  simply  allowed  to  settle  by  gravity.  The  settled  solutions  were 
decanted,  the  supernatant  was  assayed  for  viable  microbe  counts.  The  number  of 
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microbes  in  the  supematent,  relative  to  the  number  originally  added  to  the  granular  or 
colloidal  particle  solution,  gave  an  indication  of  the  number  of  viruses  either  inactivated 
or  adsorbed  to  the  particles.  However,  this  assay  could  not  distinguish  whether  the 
microbes  were  killed  or  simply  adsorbed  to  the  solids.  Therefore,  a 5 % beef  extract 
solution  was  added  to  the  solids  in  an  attempt  to  elute  any  adsorbed  microbes.  The  beef 
extract-particle  mixtures  were  agitated  for  10  minutes,  and  the  supernatant  was  again 
assayed  for  viable  microbes  to  determine  if  microorganisms  may  have  been  adsorbed  but 
not  inactivated.  In  addition,  a mixture  of  deionized  water  and  AEM  5700  coated  silica 
solution,  which  was  free  of  microorganisms,  was  agitated,  and  the  supernatant  was  tested 
for  its  antimicrobial  activity  in  a manner  analogous  to  that  described  above. 

Sand  that  was  coated  with  metal  hydroxides  and  AEM  5700,  silica  that  was 
coated  with  AEM  5700  and  diatomaceous  earth  that  was  coated  with  metal  hydroxides  all 
showed  significant  increase  of  MS-2  inactivation  compared  to  controls,  as  seen  in  Figure 
6-1 . However,  Figure  6-2  shows  that  AEM  5700  and  metal  hydroxide  coated  sand  were 
not  as  effective  in  improving  the  inactivation  of  E.  coli,  but  smaller  particles 
(diatomaceous  earth  and  silica)  were  able  to  significantly  improve  the  inactivation  of  E. 
coli  coated  with  either  metal  hydroxides  or  AEM  5700,  respectively.  Curiously,  coated 
and  uncoated  spherical  glass  particles  inactivated  less  MS-2  and  E.  coli  than  virtually  any 
other  media.  The  effect  of  substrate  in  creating  a suitable  antimicrobial  surface  seems  to 
be  significant  but  not  well  understood.  Tables  6-1  and  6-2  illustrate  that  beef  extract 
solutions  were  unable  to  elute  significant  numbers  of  MS-2  or  E.  coli  from  metal 
hydroxide  or  siloxane  modified  media,  indicating  that  microbes  were  likely  inactivated 
instead  of  simply  adsorbed.  The  AEM  5700  supernatant  control  showed  no  antimicrobial 
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Figure  6-1 . Batch  inactivation  of  MS-2  by  colloidal  and  granular  particles. 
Metal  hydroxide  and  AEM  5700  coated  materials  outperform  uncoated 
media  except  in  the  case  of  coated  glass  beads.  The  last  two  items  represent 
controls  of  unbound  AEM  5700  added  to  MS-2  and  rinsate  of  AEM  5700 
coated  particles,  respectively. 


activity,  indicating  that  the  siloxane  was  not  leached  from  the  surface,  which  was  in 
keeping  with  previous  studies  (Clarkson  and  Evans,  1993;  Isquith  et  al,  1972).  Also,  the 
unbound  AEM  5700  that  was  combined  with  the  microbes  reveal  that  the  quaternary 
ammonium  siloxane  itself  has  strong  antimicrobial  activity.  The  batch  inactivation 
studies  suggest  that  both  the  quaternary  ammonium  siloxane  AEM  5700  and  the  mixture 
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of  ferric  and  aluminum  hydroxide  are  promising  as  coatings  for  promoting  the  filtration 
of  bacteria  and  viruses. 


! 
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Figure  6-2.  Batch  inactivation  of  E.  coli  by  colloidal  and  granular 
particles.  Metal  hydroxide  coated  diatemaceous  earth  and  AEM  5700 
coated  silica  and  sand  outperformed  uncoated  media.  The  last  two  items 
represent  controls  of  unbound  AEM  5700  added  to  E.coli  and  rinsate  of 
AEM  5700  coated  particles,  respectively. 
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Table  6-1 . Inactivation  and  subsequent  elution  of  E.  coli  in  batch  tests. 


Media 

Initial  Concentration 
(cfu/ml) 

Supernatant 

Concentration 

Eluant 

Concentration 

AEM  Coated  Sand 

3.6x10* 

2.9x1 04 

2.0x1 04 

Control  Sand 

3.6x1 06 

7.2x1 06  1 

None  detected 

AEM  Coated  Silica 

3.6x1 06 

None  detected 

None  detected 

Control  Silica 

3.6x1 06 

2.7x10* 

9.0x10* 

Metal  hydroxide 
diatomaceous  earth 

6.1xl05 

2.4x1 02 

- 

No  particles* 

4.8x10' 

5.1x10' 

- 

AEM  only** 

5.6x10* 

- 

- 

* This  represents  the  control  group  in  which  only  the  supematnat  from  coated  silica 
solution  was  charged  with  E.  Coli. 

**  This  control  was  for  E.  coli  added  to  free  AEM  5700. 


Table  6-2.  Inactivation  and  subsequent  elution  of  MS-2  in  batch  tests. 


Media 

Initial  Concentration 
(pfu/ml) 

Supernatant 

Concentration 

Eluant 

Concentration 

Coated  Sand 

2.9x10* 

None  detected 

None  detected 

Control  Sand 

2.9x10* 

2.5x10* 

2.2x10* 

Coated  Silica 

2.9x10* 

8.2x10* 

3.0x10* 

Control  Silica 

2.9x10* 

2.2x10* 

2.9x1 04 

Metal  hydroxide 
diatomaceous  earth 

7.1xl07 

1.2xl03 

- 

No  particles* 

7.1x10' 

9.5x10' 

- 

AEM  only** 

2.9x10* 

- 

- 

* This  represents  the  control  group  in  which  only  the  supematnat  from  coated  silica 
solution  was  charged  with  MS-2. 

**  This  control  was  for  E.  coli  added  to  free  AEM  5700. 


6-3  Column  Filtration  of  Viruses 

Four  identical  filtration  columns,  as  described  in  5-1,  were  filled  with  either 
Ottawa  sand,  AEM  5700  coated  Ottawa  sand,  activated  carbon  or  metal  hydroxide  coated 
Ottawa  sand.  Viral  solutions  were  passed  top-to-bottom  through  the  columns  at  a 
superficial  velocity  of  about  1 mm/s,  which  is  typical  of  industrial  and  other  water 


105 


treatment  processes.  The  same  stock  solution  was  fed  to  all  four  columns.  Pressure 
drops  were  measured  in  each  of  the  four  columns  using  a mercury  manometer.  Samples 
were  taken  for  assay  from  the  inlet  and  from  the  exit  of  each  column  periodically  for 
about  two  hours  or  about  70  pore  volumes,  ensuring  steady-state  conditions.  Virus 
removals  for  each  column  and  organism  were  determined  by  viral  assays,  and 
corresponding  filter  coefficients  are  reported  in  Figure  6-3.  Metal  hydroxide  coated  sand 
was  found  to  remove  52  + 7%  and  46  + 9%  of  MS-2  and  PRD-1,  respectively,  whereas 
unmodified  sand  columns  were  seen  to  remove  -3  + 4%  and  6 + 9%,  respectively.  A 
two-factor  analysis  of  variance  revealed  a significant  difference  (a  = 0.05)  in  percent 
removals  between  different  sand  columns.  Tukey's  test  was  performed  with  a = 0.05  to 
directly  compare  the  percent  removals  of  pairs  of  columns.  Metal  hydroxide  coated 


Uncoated  Unmodified  AEM  5700  Metal 

Sand  Carbon  Coated  Sand  Hydroxide 

Sand 


Figure  6-3.  Column  filtration  of  bacteriaphages  MS-2  and  PRD- 1 . Iron  and 
aluminum  hydroxide  coated  sand  outperforms  all  other  media 
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Ottawa  sand  removed  more  MS-2  and  PRD-1  than  both  activated  carbon  and  uncoated 
Ottawa  sand  and  more  PRD-1  than  activated  carbon  and  AEM5700  coated  Ottawa  sand. 
AEM5700  coated  sand  removed  significantly  more  MS-2  than  unmodified  sand  but  not 
significantly  more  PRD-1,  nor  did  AEM  5700  coated  columns  remove  more  viruses  than 
activated  carbon  columns.  The  mediocre  performance  of  AEM  5700  coated  filter  media 
is  in  contrast  to  its  excellent  performance  in  batch  tests.  One  possible  explanation  for  this 
discrepancy  may  be  that  the  decline  in  positive  zeta  potential  results  in  an  electrostatic 
barrier  between  viruses  and  coated  media.  Another  explanation  is  that  the  hydrophobic 
nature  of  the  siloxane-coated  media  may  affect  the  transport  of  the  solution,  perhaps 
resulting  in  fluid  entrainment  in  the  column. 

6-4  Economic  Feasibility  of  Metal  Hydroxide  Coated  Filter  Media 

The  chemical  costs  of  coating  filter  media  with  iron  aluminum  hydroxides  were 
estimated  based  on  observed  chemical  requirements  for  coating  media  in  the  laboratory, 
together  with  costs  of  chemicals  in  bulk  quantities  (Table  6-3).  The  total  coating  cost  for 
sand  was  7.0  cents  per  kg,  whereas  the  coating  cost  for  activated  carbon  was  55  cents  per 
kg.  These  costs  were  used  in  determining  the  added  price  to  the  consumer  of  water 
treated  by  coated  filter  media. 

The  lower  of  the  mean  percent  removals  (46%  achieved  with  PRD  1 ) was  chosen 
to  be  representative  of  virus  removals  obtained  by  coated  sand.  Given  this  level  of 
performance  in  a 0.5  m deep  column,  a corresponding  removal  of  70%  may  be  expected 
in  a 1 .0  m deep  column  based  on  first-order  removal  (Tien,  1989).  The  cost  of  coating 
filter  media  to  achieve  this  removal  can  be  based  on  a unit  cross-sectional  area  of  1 .0  m2 
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(normal  to  the  downward  flow  direction  of  the  water  being  treated).  The  chemical  cost  to 
coat  the  media  contained  in  this  unit  area  would  be  $120.57  (Table  6-4).  The  volume  of 
water  treated  in  1 month  by  this  unit  area,  based  on  a superficial  velocity  of  1 .4  mm/s  [ 2 
gal/(fr  min)  ] and  a 5%  loss  for  backwashing,  would  be  3498  m3.  The  unit  chemical 
cost  for  water  treatment  with  coated  media  is  thus  computed  to  be  3.4  cents  per  cubic 
meter.  This  compares  to  a typical  delivered  water  price  in  U.S.  municipal  systems  of  40 
cents  per  cubic  meter. 

The  cost  of  water  treatment  by  a conventional  household  faucet  filter  as  compared 
to  the  chemical  cost  of  coating  the  activated  carbon  in  such  a filter  is  determined  in  Table 
6-5.  Based  on  a retail  cost  of  $10  per  filter,  a lifetime  of  1 month,  and  a daily  throughput 
of  8 L,  the  unit  cost  of  water  treated  by  a faucet  filter  was  estimated  to  be  4.1  cents  per 
liter.  Coating  the  activated  carbon  contained  in  one  faucet  filter  would  require  2.9  cents 
worth  of  chemicals,  amounting  to  an  additional  cost  of  0.01  cents  per  liter. 


Table  6-3.  Chemical  costs  of  coating  sand  and  activated  carbon 


FeCl36H20 

AICI36H2O 

NH4OHa 

Total 

Bulk  chemical  costb 

$0. 33/kg 

$ 1.76/kg 

$0.25/kg 

— 

Amount  used  to 
coat  1 kg  media 

Sand 

16.9  g 

30.2  g 

45  g 

— 

Activated 

carbon 

135  g 

231  g 

346  g 

— 

Cost  in  cents 
per  kg  media 

Sand 

0.5 

5.3 

1.2 

7.0 

Activated 

carbon 

4.4 

41 

9.5 

55 

a28-30%  as  NH3  with  a specific  gravity  of  approximately  0.9 
bChemical  Market  Reporter  (200 1 ) 
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Table  6-4.  Cost  of  using  coated  sand  for  municipal/industrial  filtration 


Filter  coefficient 

1.2  m'1 

Bed  depth 

1.0  m 

Initial  virus  removal 

70% 

Porosity 

0.35 

Bed  cross-sectional  area 

1.0  m2 

Bed  depth 

1.0  m 

Bed  volume 

1.0  m3 

Media  weight* 

1722.5  kg 

Unit  coating  cost 

$0.07/kg 

Coating  cost  for  1 .0  m3  bed 

$120.57 

Superficial  velocity 

1 .4  mm/s 

Volume  treated  in  1 month 

3498  m3 

Unit  cost 

3.4  cents/m3 

Typical  U.S.  delivered 
municipal  water  cost 
(conventional  filter  media) 

40  cents/m3 

* Specific  gravity  of  sand  = 2.65 


The  estimated  chemical  costs  for  coating  sand  and  activated  carbon,  when  expressed 
relative  to  the  volume  of  water  treated,  are  low  relative  to  typical  consumer  costs  for 
water.  The  chemical  costs  of  coating  would  add  approximately  10%  to  the  cost  of  water 
delivered  by  municipal  systems,  whereas  coating  chemical  costs  would  add  less  than  1 % 
to  the  cost  of  water  treated  by  faucet  filters.  It  is  thus  clear  that  the  improvement  in  virus 
removal  performance  gained  by  use  of  coated  filter  media  provides  a significant  benefit 
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to  the  consumer  in  terms  of  increased  microbiological  quality  at  a negligible-to-modest 
increase  in  cost. 


Table  6-5.  Cost  of  using  coated  activated  carbon  in  faucet  filters 


Media  weight  per  filter 

52  g 

Initial  virus  removal 

98% 

Unit  coating  cost 

$0.55/kg 

Coating  cost  for  1 filter 

2.9  cents 

Daily  water  consumption 

8 L 

Volume  treated  in  1 month 

243.5  L 

Typical  retail  cost  of  conventional 
faucet  filter 

$10 

Water  cost  for  conventional 
faucet  filter  (assuming  change-out 
after  1 month) 

4.1  cents/L 

Additional  chemical  cost  for 
using  coated  media  in  faucet  filter 

0.01  cents/L 
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